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PART-A 



EXPERIMENT NO-1 



AIM: VERIFICATION OF SAMPLING THEOREM 



Sampling: Is the process of converting a continuous time signal into a discrete time 
signal. It is the first step in conversion from analog signal to digital signal. 

Sampling theorem: Sampling theorem states that "Exact reconstruction of a continuous 
time base-band signal from its samples is possible, if the signal is band-limited and the 
sampling frequency is greater than twice the signal bandwidth", 
i.e. f s > 2W, where W is the signal bandwidth. 

Nyquist Rate Sampling: The Nyquist rate is the minimum sampling rate required to 
avoid aliasing, equal to the highest modulating frequency contained within the signal. 
In other words, Nyquist rate is equal to two sided bandwidth of the signal (Upper and 
lower sidebands) 
To avoid aliasing, the sampling rate must exceed the Nyquist rate. i.e. f s > f^ 



Program: 

% Nyquist Rate Sampling. 

clc; 

fs = 1400; 

t = 0:l/fs:13/fs; 

x = cos(2*pi*400*t)+cos(2 

xm = abs (fft(x)); 

disp('xm'); 

disp (xm); 

k = 0:length(xm)-l; 

subplot (2,2,1); 

stem (100*k, xm); 

xlabel ('Hz'); 

ylabel ('Magnitude'); 

title ( f NR sampling 1 ); 

%UNDER Sampling. 

fs=1000; 

t = 0:l/fs:9/fs; 

x = cos(2*pi*400*t)+cos(2 

xml = abs(fft(x)); 

disp('xmr); 

disp (xml); 

k = 0:length(xml)-l; 

subplot(2,2,2); 

stem(100*k,xml); 

xlabel( f Hz f ); 

ylabel( f Magnitude 1 ); 



% Clear screen 

% Sampling frequency, fs = 1400Hz 
% Number of samples 
*pi*700*t); % Input signal 

% Determine the absolute FFT of the input signal 

% Displays xm on command window 

% Displays values of xm on command window 

% Number of samples to be plot on x-axis 

% Divide the figure window to plot the output 

% Plot the output 

% Name x-axis as "Hz" 

% Name y-axis as "Magnitude" 

% Title is "NR Sampling" 



% Sampling frequency, fs = 1000Hz 
% Number of samples 
*pi*700*t); % Input signal 

% Determine the absolute FFT of the input signal 

% Displays xml on command window 

% Displays values of xml on command window 

% Number of samples to be plot on x-axis 

% Divide the figure window to plot the output 

% Plot the output 

% Name x-axis as "Hz" 

% Name y-axis as "Magnitude" 



title( f UNDER sampling'); % Title is "UNDER Sampling" 

%OVER Sampling. 

fs = 2000; % Sampling frequency, fs = 2000Hz 

t = 0: l/fs:20/fs; % Number of samples 

x = cos(2*pi*400*t)+cos(2*pi*700*t); % Input signal 

xm2 = abs(fft(x)); % Determine the absolute FFT of the input signal 

disp('xm2'); % Displays xm2 on command window 

disp (xm2); % Displays values of xm2 on command window 

k = 0:length(xm2)-l ; % Number of samples to be plot on x-axis 

subplot(2,2,3); % Divide the figure window to plot the output 

stem(100*k,xm2); % Plot the output 

xlabel( f Hz f ); % Name x-axis as "Hz" 

ylabel( f Magnitude 1 ); % Name y-axis as "Magnitude" 

title( f OVER sampling'); % Title is "OVER Sampling" 

Output : 

xm 
Columns 1 through 7 
0.0000 0.0000 0.0000 0.0000 7.0000 0.0000 0.0000 

Columns 8 through 14 
14.0000 0.0000 0.0000 7.0000 0.0000 0.0000 0.0000 

xml 
Columns 1 through 7 
0.0000 0.0000 0.0000 5.0000 5.0000 0.0000 5.0000 

Columns 8 through 10 
5.0000 0.0000 0.0000 

xm2 
Columns 1 through 7 
2.0000 2.0741 2.3496 3.1607 11.5127 0.4052 1.8998 

Columns 8 through 14 
8.6165 4.2355 1.2552 0.3357 0.3357 1.2552 4.2355 

Columns 15 through 21 
8.6165 1.8998 0.4052 11.5127 3.1607 2.3496 2.0741 



Graph: 




EXPERIMENT NO-2 

AIM: IMPULSE RESPONSE OF A GIVEN SYSTEM 

A discrete time system performs an operation on an input signal based on a predefined 
criteria to produce a modified output signal. The input signal x(n) is the system 
excitation, and y(n) is the system response. The transform operation is shown as, 



x(n) 




y(n) = T[(x(n)] 

► 



If the input to the system is unit impulse i.e. x(n) = 5(n) then the output of the system is 
known as impulse response denoted by h(n) where, 

h(n) = T[5(n)] 

we know that any arbitrary sequence x(n) can be represented as a weighted sum of 
discrete impulses. Now the system response is given by, 

00 



y(n) = T[x(n)] = T 
For linear system (1) reduces to 



2>(k) 5(n-k)] 



..(1) 



K=-00 



00 



y(n) = £x(k) T[8(n-k)] 
k=-oo 



...(2) 



The response to the shifted impulse sequence can be denoted by h(n, k) is denoted by, 
h(n, k) = T[S(n-k)] ...(3) 



For a time-invariant system 

h(n, k) = h(n-k) 

Using (4) in (3) we obtain, 

T[5(n-k)] = h(n-k) 

Using (5) in (2) we get, 



...(4) 



•(5) 



GO 



y(n) = £x(k) h(n . k) 

k=-oo 



...(6) 



For a linear time-invariant system if the input sequence is x(n) and impulse response 

h(n) is given, we can fine output y(n) by using eqn (6). 

Which is known as convolution sum and can be represented by y(n) = x(n) * h(n) 

For Example let's find out an impulse response of a difference equation. 

The general form of difference equation is, 

M M 

y(n) = £a k y(n . k) + £b k x(n . k) . . . (7) 

k=l k=0 

For input x(n) = 5(n) 

M 

y(n) = Ypk x(n-k) = for n > M 

k=0 
Now (7) can be written as, 

N 

y(n) = £a k y(n-k) = a<>=l ...(8) 

k=0 
The solution of eqn (8) is known as homogeneous solution. The particular solution is 
zero since x(n) = for n > 0, that is 

y P (n) = o 

Therefore we can obtain the impulse response by solving the homogeneous equation and 
imposing the initial conditions to determine the arbitrary constants. 

Example: 

Let's take a second order difference equation 

y(n) - (1/6) y(n-l) -(1/6) y(n-2) = x(n) . . .(9) 

For impulse response the particular solution y p (n) = 

So, y(n) = y h (n) ...(10) 

Lety h (n) = r ...(11) 

Substituting (1 1) in (9) we get, 

r - ( 1 16) X nA - ( 1 16) r~ 2 = ( x(n) = for homogeneous solution) 

or 

X 2 -{\I6)X- (1/6) = ...(12) 



The roots of the characteristic equation are, 
A,i= 1/2,^ = -1/3 

The solution y h (n) = d (l/2) n + C 2 (-l/3) n . . .(13) 

For impulse x(n) = 8(n); x(n) = for n > and x(0) = 1 

Substituting the above relations in eqn (9) we have, 

For n = 0, 

y(0) - (1/6) y(-l) -(1/6) y(-2) = x(0) = 1 

i.e. y(0) = l. ...(14) 

For n = 1 , 

y(l)-(l/6)y(0)-(l/6)y(-l) = x(l) 

y(l)-(l/6) = 

y(l)=l/6. ...(15) 

From eqn (13) 

y(0) = C 1 + C 2 ...(16) 

y(l) = (l/2)C 1 -(l/3)C 2 ...(17) 

Comparing equations (14), (15), (16) and (17) we get, 
Ci + C 2 =l 
(l/2)Ci-(l/3)C 2 =(l/6) 

Solving for Ci and C2 we get, 
Ci = 3/5, C 2 = 2/5 

Therefore, the solution 

y(n) = (3/5)(l/2) n + (2/5)(-l/3) n ...(18) 

Now, let's find out impulse response. 
We know that, 

y(n) = x(n) * h(n) 

h(n) = y(n) / x(n) ...(19) 



Compute y(n) for various values of n 
Substitute n = 0, 1, 2, 3 in eqn (18) we get, 

y(0)= 1.0000 
y(l) = 0.1667 
y(2) = 0.1944 
y(3) = 0.0602 

So, y(n) = {1.0000, 0.1667, 0.1944, 0.0602} 
We know that, x(n) = 1 

So, the impulse response for the given difference equation using eqn (19) is 
h(n) = {1.0000, 0.1667, 0.1944, 0.0602} 

Program: 

clc; % Clear screen 

nr = [ 1 ] ; % Numeartor co-efficients 

dr = [ 1 ,-1/6,-1/6] ; % Denominator co-efficients 

h = impz(nr,dr,4); % Computes the impulse response of a given 

% difference equation, returns the co-efficients 
disp( f Impulse response h(n): f );% Displays impulse response 
disp(h); % Displays impulse response on command window 

subplot(2,l,l); % Divide the figure window to plot the output 

stem(h); % Displays the impulse response 

title( f Impulse response h(n): ? );% Title as impulse response 
x = [ 1 ] ; % Input co-efficients 

y = conv (h, x); % Convolution of input and impulse co-efficients to get the 

% output 
disp( f Output y(n): f ) % Displays Output y(n) 

disp(y); % Displays the output on command window 

subplot(2,l,2); % Divide the figure window to plot the output 

stem(y); % Plots the output 

title( f Output y(n) for given x(n) f ); % Title as "Output y(n) for given x(n)" 

Output: 

Impulse response: 
1.0000 
0.1667 
0.1944 
0.0602 
Output y(n): 

1.0000 

0.1667 

0.1944 

0.0602 
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Graph: 
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EXPERIMENT NO-3 

AIM: TO IMPLEMENT LINEAR CONVOLUTION OF TWO GIVEN 
SEQUENCES 

Theory: Convolution is an integral concatenation of two signals. It has many 
applications in numerous areas of signal processing. The most popular application is the 
determination of the output signal of a linear time-invariant system by convolving the 
input signal with the impulse response of the system. 

Note that convolving two signals is equivalent to multiplying the Fourier Transform of 
the two signals. 

Mathematic Formula: 



The linear convolution of two continuous time signals x(t) and h(t) is defined by 

i 

For discrete time signals x(n) and h(n), is defined by 



1=: V 



?<p) = i(h) a m = 2^ j$) JKk - £) 



Where x(n) is the input signal and h(n) is the impulse response of the system. 

In linear convolution length of output sequence is, 
length(y(n)) = length(x(n)) + length(h(n)) - 1 

Graphical Interpretation: 

• Reflection of h(k) resulting in h(-k) 

• Shifting of h(-k) resulting in h(n-k) 

• Element wise multiplication of the sequences x(k) and h(n-k) 

• Summation of the product sequence x(k) h(n-k) resulting in the convolution 
value for y(n) 

Example: 

x(n)={l,2,3, 1} 
h(n)= {1,1,1} 

length(y(n)) = length(x(n)) + length(y(n)) - 1 
=4+3-1=6 



12 



x(k) 



h(k) 



-4-3-2-101234567 

n=0; h(-k ) 

1 

-4-3-2-101234567 

n=l;h(l-k) 
1 

-4-3-2-101234567 

n=2; h(2-k) 

1 

-4-3-2-101234567 



n=3; h(3-k) 



.4-3-2-101234567 



n=4; h(4-k) 



-4-3-2-101234567 



n=5; h(5-k) 



-4-3-2-101234567 



1 1 1 





oo 



y(0) = X x(k) h(-k) = 1 
k=-oo 



00 



y(l) = 2>(k)h(l-k) = 1+2=3 

k=- oo 



00 



y(2) =X x(k) h(2-k) = 1+2+3=6 
k=-oo 



00 



y(3) =X x(k) h(3-k) = 2+3+1=6 
k=-oo 



GO 



y (4) = £ x (k) h(4-k) = 3+1=4 
k=-oo 



GO 



y(5) = I x(k) h(5-k) = 1 
k=-oo 
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6 6 



y(n)={l,3,6,6,4, 1} 



-4-3-2-101234567 



Program: 



clc; 

xl = input( f Enter the 1st seq: f ); 

x2 = input( f Enter the 2nd seq: f ); 

y = conv(xl, x2); 

disp(The linear convolution of two 

disp(y); 

n = 0:length(y)-l; 

stem(n, y); 

xlabel(Time f ); 

ylabel( f Magnitude f ); 

title( f Linear convolution 1 ); 



% Clear screen 
% Get the first sequence 
% Get the second sequence 
% Convolve two signals 
sequences: 1 ); % Displays the result 
% Displays the result on command window 
% Defines the length for x-axis 
% Plots the output 
% Name the x-axis as Time 
% Name the y-axis as Magnitude 
% Title as "Linear convolution" 



Output: 

Enter the 1st seq:[l 2 3 1] 
Enter the 2nd seq:[l 1 1] 
The linear convolution of two sequences: 
13 6 6 4 1 
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Graph: 
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EXPERIMENT NO-4 



AIM: TO IMPLEMENT CIRCULAR CONVOLUTION OF TWO GIVEN 



SEQUENCES 



Circular convolution: 



Let Xi(n) and X2(n) are finite duration sequences both of length N with DFT's Xi(k) 
and X2(k). Convolution of two given sequences Xi(n) and X2(n) is given by the 
equation, 

x 3 (n) = IDFT[X 3 (k)] 
X 3 (k) = Xi(k) X 2 (k) 

N-l 

x 3 (n) = X x i( m ) x 2 ((n-m)) N 

m=0 

Example: 

Let's take Xi(n) = {1, 1, 2, 1} and 
x 2 (n)={l,2,3,4} 

Arrange Xi(n) and X2(n) in circular fashion as shown below. 



xi(m) 



x 2 (m) 




xi(3) 
1 




x 2 (3) 
4 
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To get X2(-m), rotate X2(m) by 4 samples in clockwise direction. 
x 2 (-m) 




x 2 (l) 

2 

x 3 (0)= x 1 (m)x 2 (-m) 

= x,(0) x 2 (0) + x,(l) x 2 (3) + x,(2) x 2 (2) + x,(3) x 2 (l) 

= 1+4 + 6+2 
x 3 (0) = 13 

Keep Xi(m) constant and rotate x 2 (-m) once to compute further values. 

To get X3(l) rotate x 2 (-m) by one sample in anti-clockwise direction 

x 2 (l-m) 




x 2 (2) 
3 

x 3 (l) = x 1 (m)x 2 (l-m) 

= x,(0) x 2 (l) + x,(l) x 2 (0) + x,(2) x 2 (3) + x,(3) x 2 (2) 

=2+1+8+3 
x 3 (l) = 14 
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To get X3(2) rotate X2(l-m) by one sample in anti-clockwise direction 
x 2 (2-m) 



x 2 (0) 




x 3 (2) = x^m) x 2 (2-m) 

= x,(0) x 2 (2) + xi(l) x 2 (l) + xi(2) x 2 (0) + xi(3) x 2 (3) 

= 3+2 + 2+4 
x 3 (2) = 1 1 

To get x 3 (3) rotate x 2 (2-m) by one sample in anti-clockwise direction 

x 2 (3-m) 




x 2 (0) 
1 

x 3 (3) = xi(m) x 2 (3-m) 

= x,(0) x 2 (3) + x,(l) x 2 (2) + xi(2) x 2 (l) + x,(3) x 2 (0) 

=4+3+4+1 

x 3 (3)=12 

The convoluted signal is, 
x 3 (n)={13, 14, 11,12} 



% Clear screen 
% Get the first sequence 
% Get the second sequence 
% Returns the length of first sequence 
% Returns the length of second sequence 
% Get the Maximum length out of two signals 



% Pad zeros to the first sequence if N2>N1 
% Pad zeros to the second sequence is N1>N2 



Program: 

clc; 

xl = input( f Enter the first sequence:'); 
x2 = input( f Enter the second sequence: 1 ); 
Nl=length(xl); 
N2 = length(x2); 
N = max(Nl,N2); 
N3=N1-N2; 

% Get equal length sequence 
if(N3>0) 

xl = [xl,zeros(l,N3)]; 
else 

x2 = [x2,zeros(l,-N3)]; 
end 

% Computation of circular convolution 
form=l:N 
y(m) = 0; 

forn=l:N 
i=m-n+l; 
if(i<=0) 

i=N+i; 
end 
y(m) = y(m)+xl(n)*h(i); 

end 
end 

disp(The circular convolution of two given sequences is: f ); 

disp(y); % Displays the result on command window 

n=0:N- 1 ; % Get the length for x-axis 

stem(n, y); % Plot the result 

xlabel(Time f ); % Name the x-axis as "Time" 

ylabel( f Magnitude 1 ); % Name the y-axis as "Magnitude" 

title( f Circular convolution 1 ); % Title as Circular convolution 

Output: 

Enter the first sequence: [112 1] 
Enter the second sequence: [12 3 4] 
The circular convolution of two given sequences is: 
13 14 11 12 
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Graph: 
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EXPERIMENT NO-5 

AIM: AUTOCORRELATION OF A GIVEN SEQUENCE AND VERIFICATION 
OF ITS PROPERTIES 

Correlation: Correlation determines the degree of similarity between two signals. If the 
signals are identical, then the correlation coefficient is 1 ; if they are totally different, the 
correlation coefficient is 0, and if they are identical except that the phase is shifted by 
exactly 180°(i.e. mirrored), then the correlation coefficient is -1. 

Autocorrelation: The Autocorrelation of a sequence is correlation of a sequence with 
itself. The autocorrelation of a sequence x(n) is defined by, 



R xx (k) = X x(n) x(n-k) k = 0, ± 1 , ±2, ±3 

n=-oo 
Where k is the shift parameter. 
Or equivalently 



Rxx (k) = X x(n+k) x(n) k = 0, ±1 , ±2, ±3 

n=-oo 

Example: 

Let's take a sequence x(n) = {1, 2, 3, 4} 



x(n) 



-4-3-2-101234567 



k=0; x(n) 



00 



Rxx(O) = I x(n) x(n) = 4 



n=-oo 



-4-3-2-101234567 
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k=l;x(n+l) 



00 



R X x(l) = Zx(n)x(n+l) = 8+3=11 



n=- co 



-4-3-2-101234567 



k=2; x(n+2) 2 



co 



Rxx(2) =X x(n) x(n+2) = 12+6+2=20 



n=-co 



-4-3-2-101234567 



k=3; x(n+3) 2 
1 



co 



R xx (3) =2 x(n) x(n+3) = 16+9+4+1=30 



n=-oo 



-4-3-2-101234567 



k=4; x(n+4) 2 
1 



-4 -3 -2 - 



co 



R xx (4) = X x ( n ) x ( n+4 ) = 12+6+2=20 



n=-co 



12 3 4 5 6 7 



co 



k=5;x(n+5) _ 4 

3 
2 
1 



R xx (5) = Zx(n)x(n+5) = 8+3=11 



n=-oo 



-4-3-2-101234567 



k=6;x(n+6) 4 

3 
2 
1 



-4-3-2-101234567 



oo 



R xx (6) = I x(n) x(n+6) = 4 



n=-co 
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30 



20 



11 



20 



11 



-4-3-2-101234567 



R xx (k)={4, 11,20,30,20, 11,4} 



Program: 

clc; 

x = input( f Enter the input sequence:'); 

Rxx = xcorr(x); 

disp( f Auto correlated sequence, Rxx: f ); 

disp(Rxx); 

t = 0:length(Rxx)-l; 

stem(t,Rxx); 

xlabel(Time f ); 

ylabel( f Magnitude 1 ) ; 

title( f Auto-correlation of the given sequence: 



Clear screen 

Get the input sequence 

Returns auto-correlated sequence 

Display the result 

Display the result on command window 

Length to plot the graph 

plots the result on figure window 

xlabel 

ylabel 

% Title 



Output: 

Enter the input sequence: [12 3 4] 
Auto correlated sequence, Rxx: 

4.0000 11.0000 20.0000 30.0000 



20.0000 11.0000 4.0000 
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Graph: 




Properties of Autocorrelation: 

i. Periodicity: R xx (ko) = R X x(0) then R xx (k) is periodic with period ko 

2. Autocorrelation function is symmetric, i.e. R xx (m) = R xx (-m) 

3. Mean square value: autocorrelation function at k=0, is equal to mean square 
value of the process. Rxx(O) = E{|x(n)| }>0 



24 



EXPERIMENT NO-6 

AIM: CROSS-CORRELATION OF A GIVEN SEQUENCE AND VERIFICATION 
OF ITS PROPERTIES 

Correlation: Correlation determines the degree of similarity between two signals. If the 
signals are identical, then the correlation coefficient is 1 ; if they are totally different, the 
correlation coefficient is 0, and if they are identical except that the phase is shifted by 
exactly 180°(i.e. mirrored), then the correlation coefficient is -1. 

Cross-correlation: When two independent signals are compared, the procedure is known 
as cross-correlation. It is given by, 



R xy (k) = X *(n) Y(n-k) k = 0, ± 1 , ±2, ±3 

n=-oo 
Where k is the shift parameter. 
Or equivalently 



Ryx(k) = Zy( n+k ) x ( n ) 



n=-oo 



Comparing above two equations, we find that, 

R xy (k) = R yx (-k) 

Where R yx (-k) is the folded version of R xy (k) about k = 0. 

So, we can write Cross correlation of the two sequences is given by, 

00 

Rx y (k) = Ix(n)y[-(k-n)] 

n=-oo 

R xy (k) = x(k)*y(-k) ...(1) 

Equation (1) shows that cross correlation is the essentially the convolution of two 
sequences in which one of the sequences has been reversed. 



25 



Example: 

Let's take x(n) ={1,1,1,1} and y(n) ={1,2, 3, 4} 
x(n) 

1111 



.4-3-2-101234567 



k=0; y(-n) 



-4-3-2-101234567 



k=l;y(l-n) 



-4-3-2-101234567 



k=2; y(2-n) 



-4-3-2-101234567 



k=3;y(3-n) 



00 



Rx y (0) = I x(n) y(-n) = 4 



n=-oo 



00 



R X y(l) = Zx(n)y(l-n) = 4+3=7 



n=- oo 



00 



R xy (2) =1 x(n) y(2-n) = 4+3+2=9 



n=-oo 



00 



R xy (3) =2 x(n) y(3-n) = 4+3+2+1=10 



n=-oo 



-4-3-2-101234567 



k=4; y(4-n) 
4 



-4 -3 -2 - 



oo 



R xy (4) = I x(n) y(4-n) = 3+2+1=6 



n=-oo 



12 3 4 5 6 7 
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k=5; y(5-n) 



00 



Rx y (5) = Ix(n)y(5-n) = 2+l=3 



-4-3-2-101234567 



k=6; y(6-n) 



-4-3-2-101234567 



10 



n=-oo 



oo 



R xy (6) = I x(n) y(6-n) = 1 



n=-oo 



R xy (k)= {4, 7, 9,10,6,3,1} 



-4-3-2-101234567 
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Program: 

clc; 

x = input( f Enter the first sequence: 1 ); 

y = input( f Enter the second sequence: 1 ); 

Rxy = xcorr(x,y); 

disp( f Cross correlated sequence, Rxy: f ); 

disp(Rxy); 

t = 0:length(Rxy)-l; 

stem(t, Rxy); 

xlabel(Time f ); 

ylabel( f Magnitude 1 ); 



% 



Clear screen 

Get the first sequence 

Get the second sequence 

Returns cross correlated sequence 

Display the result 

Display the result on command window 

Length to plot the graph 

plots the result on figure window 

xlabel 

ylabel 



title( f Cross correlation of the given two sequences: 1 ); % Title 

Output: 

Enter the first sequence: [1111] 
Enter the second sequence: [12 3 4] 
Cross correlated sequence, Rxy: 

4.0000 7.0000 9.0000 10.0000 6.0000 3.0000 1.0000 



Graph: 



Cross correlation of the given two sequences: 
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Properties of cross correlation: 

1 . R xy (k) is always a real valued function which may be a positive or negative. 

2. R xy (-k) = R yx (k) 

3. |R iy (k)| 2 <IU0)R y y(0) 

4. |R xy (-k)| <[l/2] [R xy (-k) + R xy (-k)] 

5. When R xy (k) = 0, x(n) and y(n) are said to be uncorrected or they said to be 
statistically independent. 

6. R xy (k) may not be necessarily have a maximum at k=0 nor R xy (k) an even 
function. 
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EXPERIMENT NO-7 
AIM: TO SOLVE A GIVEN DIFFERENCE EQUATION 

A General N order Difference equations looks like, 

y[n] + aiy[n-l] + a 2 y[n-2] + . . .+ a N y[n-N) = 

b x[n] + bix[n-l] + ....+ b M x[n-M] 

Here y[n] is "Most advanced" output sample and y[n-m] is "Least advanced" ouput 

sample. 

The difference between these two index values is the order of the difference equation. 

Here we have: n-(n-N) = N 

We can rewrite the above equation as, 

y[n] + 2>iy[n-i] = 5>ix[n-i] 

y[n] becomes, 

y[n] = -Zaiy[n-i] + £bix[n-i] 

Example: 

y[n+2] - 1.5y[n+l] + y[n] = 2x[n] 

In general we start with the "Most advanced" output sample. Here it is y[n+2]. 
So, here we need to subtract 2 from each sample argument. We get 

y[n] - 1.5y[n-l] + y[n-2] = 2x[n-2] 

«=> y[n] = 1.5y[n-l] - y[n-2] + 2x[n-2] 

Lets assume our input x[n] = u[n] = f x<0 

1 x>0 



In our example we have taken x[n] = u[n] = J 

ll 



x<0 
0<x<10 



We need N past values to solve N order difference equation. 
y[-2] = l 
y[-l] = 2 
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Compute y[n] for various values of n 
y[0] = 1.5y[-l]-y[-2] + 2x[-2] 

= 1.5*2-1+2*0 
y[0] = 2 

y[l] = 1.5y[0]-y[-l]+2x[-l] 
= 1.5*2-2 + 2*0 

y[i] = i 

y[2] = 1.5y[l]-y[0] + 2x[0] 

= 1.5*1-2 + 2*1 
y[2] = 1.5 

And so on... 



Program: 

clc; 

x=[0 0ones(l, 10)]; 
y_past = [1 2]; 
y(l) = y_past(l); 
y(2) = y_past(2); 



% Clear screen 

% Input x[n] = u[n] 

% Past values to solve difference equation 

%y[l]<=y[-2] 

%y[2]<=y[-l] 



for k = 3:(length(x)+2) % Compute y[n] for various values of n 

y(k) = 1.5*y(k-l) - y(k-2) + 2*x(k-2); 
end % End of for loop 

disp(' Solution for the given difference equation:'); 



disp(y); 

subplot( 1,1,1); 

stem (-2:(length(y)-3),y); 

xlabel('Input x[n]'); 

ylabel('Outputy[n]'); 

title(' Difference equation'); 

Output: 



% Display the result on command window 

% Divide the window to plot the result 

% Plot the result 

% Name x-axis as Input x[n] 

% Name y-axis as Output y[n] 

% Title as "Difference equation" 



Solution for the given difference equation: 
Columns 1 through 7 
1.0000 2.0000 2.0000 1.0000 1.5000 



3.2500 5.3750 



Columns 8 through 14 
6.8125 6.8438 5.4531 3.3359 



1.5508 0.9902 1.9346 
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Graph: 
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EXPERIMENT NO-8 

AIM: TO COMPUTE N-POINT DFT OF A GIVEN SEQUENCE AND TO PLOT 
MAGNITUDE AND PHASE SPECTRUM. 

Discrete Fourier Transform: The Discrete Fourier Transform is a powerful 
computation tool which allows us to evaluate the Fourier Transform X(e JC0 ) on a digital 
computer or specially designed digital hardware. Since X(e |ft) ) is continuous and 
periodic, the DFT is obtained by sampling one period of the Fourier Transform at a finite 
number of frequency points. Apart from determining the frequency content of a signal, 
DFT is used to perform linear filtering operations in the frequency domain. 

The sequence of N complex numbers Xo,..., Xjv-i is transformed into the sequence of N 
complex numbers Xo, ..., Xjv-i by the DFT according to the formula: 

N-l 

X(k) = £x(n)e- j27ink/N k = 0, 1 , . . . . N- 1 

n=0 

Example: 

Lets assume the input sequence x[n] = [ 1 10 0] 

We have, 

N-l 

X(k) = 2>(n)e- j27tnk/N k = 0, 1 , . . . . N- 1 

n=0 

For k = 0, 

3 

X(0) = £x(n) = x (0) + x(l) + x(2) + x(3) 
n=0 
X(0) = 1+1+0+0 = 2 

For k= 1, 

3 

X(l) = X x ( n ) e " J5tn/2 = x (°) + x 0) e" J5t/2 + x(2) e _j71 + x(3) e i3n/2 
n=0 

= 1 + cos(7i/2) - jsin(7i/2) 
X(l) =l-j 
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For k = 2, 
3 

X(2) = Xx(n)e" j7m = x (°) + x (!) e " j7l+ x ( 2 ) e " j271 + x ( 3 ) e " j371 
n=0 

= 1 + cos 7i-jsin7i 
X(2)= 1-1=0 

For k = 3, 
3 

X(3) = X x (n)e" j3n7l/2 = x(0) + x(l) e" j37l/2 + x(2) e" j371 + x(3) e" j97r/2 
n=0 

= 1 + cos(37i/2) - jsin(37i/2) 
X(3) =l+j 

The DFT of the given sequence is, 
X(k)={2,l-j,0,l+j} 

To find Magnitude of X(k): 

Magnitude= (a 2 +b 2 ) 1/2 

Where a and b are real and imaginary parts respectively 

To fine Phase of X(k): 

Phase=tan _1 (b/a) 



Program: 

clc; 

xl = input( f Enter the sequence:'); 

n = input( f Enter the length: 1 ); 

m = abs(fft(xl,n)); 

disp( f N-point DFT of a given sequence:'); 

disp(m); 

N = 0:l:n-1; 

subplot(2,2,l); 

stem(N,m); 

xlabel('Length'); 

ylabel('Magnitude of X(k)'); 

title('Magnitude spectrum:'); 

an = angle(fft(xl,n)); 

disp(an); 

subplot(2,2,2); 

stem(N, an); 



% Clear screen 

% Get the input sequence 

% Get the value of N 

% Computes the DFT using FFT algorithm 

% Display the results 

% Displays the result on command window 

% Decides the length to plot the results 

% Divide the figure window to plot the 

% results 

% Plots the magnitude spectrum 

% Name x-axis as "Length" 

% Name y-axis as "Magnitude of X(k)" 

% Title as "Magnitude spectrum" 

% Get the angle of the output sequence X(k) 

% Divide the figure window to plot the 

% results 

% Plots the phase spectrum 
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xlabel( f Length f ); 
ylabel(ThaseofX(k) f ); 
title(Thase spectrum: 1 ); 



% Name x-axis as "Length" 

% Name y-axis as "Phase of X(k)" 

% Title as "Phase spectrum" 



Output: 

Enter the sequence: [1 10 0] 
Enter the length: 4 
Magnitude of X(k) 

2.0000 1.4142 1.4142 



Phase of X(k) 
-0.7854 



0.7854 



Graph: 



Magnitude spectrum: 
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EXPERIMENT NO-9 

AIM: LINEAR CONVOLUTION OF TWO GIVEN SEQUENCES USING DFT 
AND IDFT 

Theory: Convolution is an integral concatenation of two signals. It has many 
applications in numerous areas of signal processing. The most popular application is the 
determination of the output signal of a linear time-invariant system by convolving the 
input signal with the impulse response of the system. 

Note that convolving two signals is equivalent to multiplying the Fourier Transform of 
the two signals. 

Mathematic Formula: 

The linear convolution of two continuous time signals x(t) and h(t) is defined by 



i 



For discrete time signals x(n) and h(n), is defined by 



1=: V 



?<p) = i(h) a m = 2^ J$) J(K - £) 



Where x(n) is the input signal and h(n) is the impulse response of the system. 

Example: 

Xl (n)={l,l,2} 
x 2 (n)={l,2} 

For linear convolution, 

Length N = Length(xi) + Length(x 2 ) - 1 

N=3+2-l=4 

Convolution of two sequences Xi(n) and X2(n) is, 
x 3 (n) = IDFT[X 3 (k)] 
x 3 (n) = IDFT[X,(k) X 2 (k)] 

Where, 

Xi(k) = DFT [ Xl (n)] 

X 2 (k) = DFT [x 2 (n)] 
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N-l 

Xi(k) = X x 1 (n)e- i2 " kD/N k= 0, 1, 2, . . . , N-l 

n=0 

Given Xi(n) = {1, 1, 2} andN=4 

3 

X 1 (0) = £x 1 (n)=l + l+2 = 4 
n=0 



X 1 (l) = Xxi(n)e- J5tn/2 = l-j-2 = -l-j 
n=0 



Xi(2) = X xi(n)e" j7in = 1-1+2 = 2 
n=0 



X 1 (3) = Xxi(n)e- j33tn/2 = l+j-2 = -l+j 
n=0 

X 1 (k)={4,-l-j,2,-l+j} 

Now, 

N-l 

X 2 (k) = X x 2 (n)e- j25tkn/N k= 0, 1 , 2, . . . , N- 1 

n=0 

Given x 2 (n) ={1,2} andN=4 

3 

X 2 (0) = X x 2 (n) = 1+2 = 3 
n=0 
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3 

X 2 (l) = Z x 2 (n) e" J7tn/2 = 1 + 2(-j) = 1 - j2 
n=0 

3 

X 2 (2) = Xx 2 (n)e- J7in =l+2(-l) = -l 
n=0 

3 

X 2 (3) = X x 2 (n) e- jM = 1 + 2(j) = 1 + j2 
n=0 

X 2 (k)={3, l-j2,-l,l+j2} 

We know that, 

X 3 (k) = Xi(k) X 2 (k) 

X 3 (k)={12,-3+j,-2,-3-j} 

Convolution of two given sequences is, 
x 3 (n) = IDFT[X 3 (k)] 

N-l 

x 3 (n) = (l/N)Xx 3 (k)^ kn/N n = 0,l,2,....,N-l 

k=0 

3 

x 3 (0) = (1/4) X X 3 (k) = (1/4) [12 - 3 +j -2 -3 -j] = 1 
k=0 

3 

x 3 (l) = (l/4)Xx 3 (k)e i7tk/2 

k=0 
x 3 (l) = (1/4) [12 + (-3 + j) j + (-2) (-1) + (-3 - j) (-j)] = 3 

3 

x 3 (2) = (1/4) X X 3 (k)d* k 
k=0 

x 3 (2) = (1/4) [12 + (-3 + j) (-1) + (-2) (1) + (-3 - j) (-1)] = 4 
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3 

x 3 (3) = (1/4) X X 3 (k)^ KkJ2 
k=0 

x 3 (3) = (1/4) [12 + (-3 + j) (-j) + (-2) (-1) + (-3 - j) (j)] = 4 



Convoluted sequence of two given sequences is, 
x 3 (n)={l,3,4,4} 



Program: 

clc; % Clear screen 

xl = input( f Enter the 1st seq: f ); % Get the first sequence 

x2 = input( f Enter the 2nd seq: ? ); % Get the second sequence 

n = length(xl) + length(x2)-l; % Get the length of the sequence 

xl = fft(xl,n); % Compute the DFT of xl using FFT algorithm 

x2 = fft(x2,n); % Compute the DFT of x2 using FFT algorithm 

y = xl .*x2; % Multiply two DFT's 

yc = ifft(y,n); % Compute IDFT using IFFT algorithm 

disp( f Linear convolution using DFT and IDFT: f ); % Display Linear convolution 

disp(yc); % Displays the result on command window 

N = 0: 1 :n-l ; % Defines the length of x-axis to plot the result 

subplot( 1,1,1); % Divide the window to plot the result 

stem(N,yc); % Plots the result 

xlabel(Time f ); % Name the x-axis as Time 

ylabel( f Magnitude f ); % Name the y-axis as Magnitude 

title( f Linear convolution using DFT and IDFT: f ); % Title 



Output: 

Enter the 1st seq: [112] 
Enter the 2nd seq: [1 2] 
Linear convolution using DFT and IDFT: 
13 4 4 
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Graph: 



Linear convolution using DFT and IDFT: 
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EXPERIMENT NO-10 

AIM: TO IMPLEMENT CIRCULAR CONVOLUTION OF TWO GIVEN 
SEQUENCES USING DFT AND IDFT 

Circular convolution: 

Let Xi(n) and X2(n) are finite duration sequences both of length N with DFT's Xi(k) 
and X2(k). Convolution of two given sequences Xi(n) and X2(n) is given by, 

x 3 (n) = IDFT[X 3 (k)] 
x 3 (n) = IDFT[X!(k) X 2 (k)] 

Where, 

Xi(k) = DFT [xi(n)] 

X 2 (k) = DFT [x 2 (n)] 

Example: 

xi(n)={l, 1,2, 1} 
x 2 (n)={l,2,3,4} 

N-l 

X!(k) = X Xl (n)e- J23lkn/N k= 0, 1, 2, ... , N-l 

n=0 

Given Xi(n) = {1, 1, 2, 1} andN=4 

3 

Xi(0) = X x i( n ) =1 + 1+2 + 1=5 
n=0 

3 

Xi(l) = X xi(n)e- jmV2 = 1 - j - 2 + j = -1 
n=0 



Xi(2) = X xi(n)e- j5m = 1-1+2-1 = 1 
n=0 



X 1 (3) = Xxi(n)e- j37tn/2 = l+j-2-j = -l 
n=0 
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Xi(k)= {5, -1,1,-1} 

Now, 

N-l 

X 2 (k) = X x 2 (n)e- j27tkn/N k= 0, 1 , 2, . . . , N- 1 

n=0 

Given x 2 (n) ={1,2, 3, 4} and N=4 

3 

X 2 (0) = X x 2 (n) =1+2 + 3+4=10 
n=0 

3 

X 2 (l) = X x 2 (n) e" J7tn/2 = 1 + 2(-j) + 3(-l) + 4(j) = -2 + j2 
n=0 

3 

X 2 (2) = X x 2 (n) e" j7tn = 1 + 2(-l) + 3(1) + 4(-l) = -2 
n=0 

3 

X 2 (3) = X x 2 (n) e- jM = 1 + 20 + 3(-l) + 4(-j) = -2 - j2 
n=0 

X 2 (k)={10,-2+j2,-2,-2-j2} 

We know that, 
X 3 (k) = X 1 (k)X 2 (k) 

X 3 (k)={50,2-j2,-2,2+j2} 

Convolution of two given sequences is, 
x 3 (n) = IDFT[X 3 (k)] 

N-l 

x 3 (n) = (1/N) X X 3 (k)d 2,tkn/N n = 0, 1, 2, . . .. , N-l 

k=0 
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x 3 (0) = (1/4) X X 3 (k) = (1/4) [50 + 2 - j2 - 2 + 2 + j2] = 13 
k=0 



jjtk/2 



x 3 (l) = (l/4) X x 3(k)^ 

k=0 
x 3 (l) = (1/4) [50 + (2 - j2) j + (-2) (-1) + (2 + j2) (-j)] = 14 



x 3 (2) = (l/4) Xx 3 (k)^ k 
k=0 

x 3 (2) = (1/4) [50 + (2 - j2) (-1) + (-2) (1) + (2 + j2) (-1)] = 11 



x 3 (3) = (1/4) X ^(k)^ 2 

k=0 
x 3 (3) = (1/4) [50 + (2 - j2) (-j) + (-2) (-1) + (2 + j2) 0)] = 12 

Convoluted sequence of two given sequences is, 
x 3 (n)={13, 14, 11,12} 



Program: 

clc; 

xl = input( f Enter 1st sequence: 1 ); 

x2 = input( f Enter 2nd sequence: 1 ) 

n = max(length(xl), length(x2)); 

xl=fft(xl,n); 

x2 = fft(x2,n); 

y = xl.*x2; 

yc = ifft(y,n); 



% Clear screen 
% Get the first sequence 
% Get the second sequence 
% Get the maximum length of the two sequences 
% Compute the DFT of the xl using FFT algorithm 
% Compute the DFT of the x2 using FFT algorithm 
% Multiply two DFT' s 

% Compute the IDFT of mutliplied sequence to get 
% the convoluted sequence 
disp( f Circular convolution using DFT and IDFT: f ); % Displays Circular convolution 
disp(yc); % Displays the result on command window 

N = 0: 1 :n-l ; % Defines the length of x-axis to plot the result 

subplot( 1,1,1); % Divide the window to plot the result 

stem(N,yc); % Plots the results 

xlabel(Time f ); % Name the x-axis as "Time" 

ylabel( f Magnitude 1 ); % Name the y-axis as "Magnitude" 

title( f Circular convolution using DFT and IDFT: f ); % Title 
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Output: 

Enter 1st sequence: [1 12 1] 
Enter 2nd sequence: [1 2 3 4] 
Circular convolution using DFT and IDFT: 
13 14 11 12 



Graph: 







Circular convolution using DFT and IDFT: 
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EXPERIMENT NO-11 

AIM: DESIGN AND IMPLEMENTATION OF FIR FILTER TO MEET GIVEN 
SPECIFICATIONS (LOW PASS FILTER USING HAMMING WINDOW) 

Finite Impulse Response (FIR) Filter: The FIR filters are of non-recursive type, 
whereby the present output sample is depending on the present input sample and previous 
input samples. 
The transfer function of a FIR causal filter is given by, 

N-l 

H(z) = £ h(n)z n 

n=0 

Where h(n) is the impulse response of the filter. 

The Fourier transform of h(n) is 

N-l 

H(d w ) = £h(n)e- jwn 

n=0 

In the design of FIR filters most commonly used approach is using windows. 

The desired frequency response H^e^) of a filter is periodic in frequency and can be 

expanded in Fourier series. The resultant series is given by, 

71 

hd(n) = (1/2tt) I H(Or dw 

-71 
And known as Fourier coefficients having infinite length. One possible way of obtaining 
FIR filter is to truncate the infinite Fourier series at n = ± [(N-l)/2] 
Where N is the length of the desired sequence. 

The Fourier coefficients of the filter are modified by multiplying the infinite impulse 
response with a finite weighing sequence w(n) called a window. 

Where w(n) = w(-n) ^ for |n| < [(N-l)/2] 
= for|n|>[(N-l)/2] 

After multiplying w(n) with hd(n), we get a finite duration sequence h(n) that satisfies 
the desired magnitude response, 

h(n) = h d (n) w(n) for |n| < [(N-l )/2] 
= for |n| > [(N-l)/2] 
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The frequency response H(e JW ) of the filter can be obtained by convolution of H^e^) 
and W(e JW ) is given by, 
71 

HCO = (1/2tt) I H^e 1 ' 9 ) W^ 9 ) d9 

-71 



Example: 

Here we design a lowpass filter using hamming window. 
Hamming window function is given by, 



w H (n) = 0.54 + 0.46 cos ((2jm)/(N-l)) 
= 

The frequency response of Hamming window is, 



for-(N-l)/2<n<(N-l)/2 
otherwise 



W H (e ]w ) = 0.54[(sin(wN/2))/(sin(w/2))] 

+ 0.23 [sin (wN/2 - ttN/N - l)/sin (w/2 - tt/N -1)] 
+ 0.23 [sin (wN/2 + tiN/N - l)/sin (w/2 + tt/N - 1)] 



Program: 

clc; 

wp = input(Tass band edge freq: f ); 
ws = input( f Stop band edge freq: f ); 
tw = ws-wp; 



% Clear screen 

% Get Passband edge frequency 

% Get Stopband edge frequency 

% Subtract PB frequency from SB frequency 



N = ceil(6.6*pi/tw)+l; 
wn = (hamming(N)); 

B = firl(N-l,wp,wn); 

disp( f Impulse response coeff= 

disp(B); 

[H,w] = freqz(B, 1,256); 



Mag = 20*loglO(abs(H)); 
plot(w/pi*pi, Mag); 
xlabel(' Frequency in radians 
ylabel( f Gain in db ? ); 



% Compute length 

% Returns N-point symmetric hamming window in 

% column vector 

% Designs (N-l) order lowpass FIR filter and 

% returns filter coefficients in length N in vector B 
f ); % Displays Impulse response coefficients 

% Displays the coefficients on command window 

% Digital filter frequency response. This function 

% returns the N-point complex frequency response 

% vector H and the N-point frequency vector w in 

% radians/sample of the filter. 

% Get the magnitude 

% Plot the Magnitude spectrum 
in terms of pi f ); % Name x-axis 

% Name y-axis as "Gain in db" 
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Output: 

Pass band edge freq: 0.05 *pi 
Stop band edge freq: 0.4*pi 
Impulse response coeff= 
Columns 1 through 7 

-0.0027 -0.0035 -0.0041 -0.0010 0.0105 0.0333 0.0666 
Columns 8 through 14 

0.1046 0.1383 0.1580 0.1580 0.1383 0.1046 0.0666 
Columns 15 through 20 

0.0333 0.0105 -0.0010 -0.0041 -0.0035 -0.0027 

Graph: 
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EXPERIMENT NO-12 

AIM: DESIGN AND IMPLEMENTATION OF IIR FILTER TO MEET GIVEN 
SPECIFICATIONS 

Basically digital filter is a linear time-invariant discrete time system. 
Infinite Impulse Response(IIR) filter : IIR filters are of recursive type, whereby the 
present output sample depends on the present input, past input samples and output 
samples. 

The impulse response h(n) for a realizable filter is, 
h(n) = for n<0 

And for stability, it must satisfy the condition, 

00 

Z i h ( n ) i < °° 

n=0 

Example: 

Let's design an analog Butterworth lowpass filter. 
Steps to design an analog Butterworth lowpass filter. 

1 . From the given specifications find the order of the filter N. 

2. Round off it to the next higher integer. 

3. Find the transfer function H(s) for Qc = lrad/sec for the value of N. 

4. calculate the value of cutoff frequency Qc 

5. find the transfer function H a (s) for the above value of Qc by substituting 
S^ (s/ Sic) in H(s). 

Program: 

clc; % Clear screen 

Pa = input(Tassband attenuation in DB: f ); % Get the passband attenuation 

Sa = input( f Stopband attenuation in DB: f ); % Get the stopband attenuation 

Fpb = input(Tassband edge frequency in Hz: f ); % Get Passband edge frequency 
Fsb = input( f Stopband edge frequency in Hz: ? ); % Get Stopband edge frequency 
fs = input( f Sampling frequency: 1 ); % Get Sampling frequency 



wp 
ws = 



■- 2*Fpb/fs; 
2*Fsb/fs; 



% Convert PB edge frequency in Hz to radians 
% Convert SB edge frequency in Hz to radians 



[N,wn] = buttord(wp,ws,Pa,Sa); 

disp( f Order: f ); 

disp(N); 

disp( f Cutoff frequency: 1 ); 

disp(wn); 



% Find cutoff frequency and order of the filter 

% Display the Order 

% Display order N on command window 

% Display Cutoff frequency 

% Display Cutoff frequency on command window 
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[b,a] = butter(N,wn); 
disp( f b= f ); disp(b); 
disp( f a- ); disp(a); 



% 



w = 0:0.01:pi; 


% 


[h,om] = freqz(b,a,w,'whole'); 


% 


m = 20*logl0(abs(h)); 


% 


an = angle(h); 


% 


subplot(2,l,l); 


% 
% 
% 


plot(om/pi,m); 


xlabel('Normalised frequency:'); 


% 


ylabel('Gain in db:'); 


% 


title(' Frequency Response:'); 


% 


subplot(2,l,2); 


% 


plot(om/pi,an); 


% 


xlabel('Normalised frequency:'); 


% 


ylabel('Phase in radians:'); 


% 


title('Phase spectrum'); 


% 


Output: 





Get Numerator and denominator coefficients of 

the filter 

Display the numerator coefficients on command 

window 

Display the denominator coefficients on 

command window 

Defines length for x-axis to plot the result 

Frequency response of the filter 

Absolute value of the frequency response vector 

Get the phase of the frequency response vector 

Divide the figure window to plot the frequency 

response 

Plot the response 

Name x-axis as "Normalised frequency" 

Name y-axis as " Gain in db" 

Title as Frequency Response 

Divide the figure window 

Plot the Phase spectrum 

Name x-axis as "Normalised frequency" 

Name y-axis as Phase in radians 

Title as Phase spectrum 



Passband attenuation in DB:4 
Stopband attenuation in DB:30 
Passband edge frequency in Hz:400 
Stopband edge frequency in Hz: 800 
Sampling frequency:2000 
Order: 
3 

Cutoff frequency: 
0.4914 



0.1600 0.4800 0.4800 0.1600 



a= 



1.0000 -0.0494 0.3340 -0.0045 
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Graph: 




-50 



T3 

- -100 
° -150 



-200 






Frequency Response 
i 1 1 \ — — 



J I I I I I I I L 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.3 0.9 

Normalised free 





Drmahsed frequency: 
Phase spectrum 



n 1 1 1 1 1 1 1 r 



0.2 0.3 0.4 0.5 0.6 0.7 O.S 0.9 
Normalised frequency: 
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PART-B 
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Manual of DSP 6713 Starter Kit 



Contents 

1 Consignment list and notes 

2 Introduction of main board 

3 The setup of USB Programmer and emulator in CCS v3.3 

4 Procedures for VTU Programs 
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Consignment list 

1 . A main board of DSP67 1 3 - INo 

2. Power Supply 5V, 3A. - INo 

3. USB Programmer with USB Cable -INo 

4. Head Phone with MIC -INo 

5. Audio Jack -3No [1+2] 

6. VTU Programs CD -INo 

7 .5 th SEM Lab DSP Manual -INo 

Introduction of main board: 

1. USB2.0 CY7C68013-56PVC, compatible with USB2.0 and USB 1.1, including 8051 

2. DSP TMS320C6713 TQFP-208 Package Device with, 4 layers board 

3. SDRAM MT48LC4M16A2 1 meg* 16 *4 bank micron 

4. FLASH AM29LV800B 8MbitlMbyte of AMD 

5. RESET chip specialize for reset with button for manually reset 

6. POWER power supply externally, special 5V, 3.3V, 1.6V chip for steady voltage 
with remaining for other devices. 

7. EEPROM 24LC64 for download of USB firmware 

8. CPLDXC95144XL 

9. AIC TLV320AIC23B sampling with 8-96KHZ, 4 channels with interface of 
headphone. 



USB 



LINE IN MIC LlNE_OUT PHONE 



'ower supply 




JTAG_XILINX 



JTAG DSP 



Board Picture: 



53 



The setup of USB Programmer and emulator in CCS v3.3 

1. Install CCS v3.3 software according to the default Custom Install. 

2. Install USB Emulator chooses to install directory from CCS v3.3 hat is if CCSv3.3 
is installed in C: / CCStudio_v3.3 directory, then install the USB emulator driver in 
this directory. 

Procedure to Setup Emulator : 

3. Open the "Setup CCStudio v3.3" 




4. Choose c67xx in the "Family". 



System Configuration 



Available Factory Boards 



WARM7 AHXD5510U5B Emulator 
WARM7 5EEDXD5510U5B Emulator 
WARM7 Sirnulatorj Big Endian 
WARM7 5imulator J Little Endian 
WARM7XD5510 Emulator 
WARM7 XD5560 Emulator 
WARM9 AHXD5510U5B Emulator 
■BARM9 5EEDXD5510U5B Emulator 



Family 



Platform 



Endianness a 



All 



C55xx 
C62xx 
C64X+ 
C64xx 



CORTEX- 
CORTEX- 



ARM9 
ARM9 



■M 



All 

AHxdsSlO., 

5EEDXds5.. 

simulator 

simulator 

xdsSlOe... 

xds560e,,, 

AHxdsSlO., 

5EEDXds5.. 



All 



big 

little 
* 
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5. Choose AHxds5 lOusb emulator in the "Platform" 





File Edit View Help 


System Configuration 


Available Factory Boards 

WC6701 Device Simulator 
WC6701 Device Simulator 
■H : C670x 5EEDXD5510U5B Emulator 
WC670xXDS510 Emulator 
WC670x XD5560 Emulator 


| Family 
||c67xx t|| 

C67xx 

C67xx 

C67xx 

C67xx 

C67xx 


Platform Endianness 

All ^||aii 


U 

i ■ 




raaEBBaniraa 

5EEDXD5510U 
simulator 
xdsSlO emulal 
xds560 emulal 


little 
big 

* 



Choose little in the "Endianness". 



File Edit View Help 


System Configuration 


Available Factory Boards 


Family 


Platform 


Endianness 


u 

: ■ 


|c67xx V||AHxds510 t| 


All _▼] 




WC671x AHXDS510USB Emulator 
WC672x AHXD5510USB Emulator 


X X 
X X 

U U 


AHxdsSlO... 
AHxdsSlO... 


All 
bjg 



7 Now you are left with two options under Available Factory Boards, Choose 
C671X AHXDS510 USB Emulator, right click and "Add to system..." 



Code Composer Studio Setup 



File Edit View Help 
System Configuration 



Available Factory Boards 



Family 



Platform Endianness 



C67xx 



AHxds5lU_^ |[AJl 



XDS510USB Emulator, 



f C672x AHXDS51UU5B Emulator^ 



flu "^Rin... * 



510... 



8. Now the Emulator and the processor both are selected under "system 
configuration". 



System Configuration 



^J My System 

B fl:C671X AH_XD55 10 Emulator 



Current Proccesor Type Driver Location 



^TM5320C671x 



C : \CC5tudio_v3 . 3\dri vers\tds6000 . d vr 
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9 Choose file and click on "save". 

10 Choose file and click on "exit" you will get a wizard as shown bellow 



Code Composer Studio Setup 




Click on yes. 



Then it launches a Code Composer Studio 
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1 1 Connect the power supply to the board (5V, 3 A): make sure that supply is there in the 
board by pressing reset button. 




1 1 After 2 minutes after power supply is turned on then connects the USB 
Programmer cable to DSP JTAG connector and Host Computer where CCS 
3.3 is installed. 
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12. Go to Debug and select the Reset Emulator 



'4/C671X AH XDS510 EmuUitor/TMS32QC671x C671x Code Composer Studit 



0M 



File Edit View Project 

■is » q | jt m 



<&T Q I S! 



0- Q GEL files 
O Projects 



m oJ3~ 



% »_; DISCONNECTED (Ur 



GEL Option Profile Tools DSP/BIOS Window Help 



Assembly/Source Stepping 


► 


Step Into 


Fll 


Step Over 


FIO 


Step Out 


Shift+Fll 



Run 
Halt 

Run Free 
Low Power Run 
Run to Cursor 
Set PC to Cursor 

: : 

Operation.. 



F5 

5hift+F5 

Alt+F5 

Ctrl+F5 

Ctrl+Alt+F5 

Ctrl+FlO 

Ctrl+Shift+FlO 

Ctrl+Shift+F5 

Ctrl+M 



Advanced Resets 


Ctrl+R 


^gj^g^^ 


Halt on Reset 


Connect 

Restore Debug State 


Alt+C 



Thread Level Debugging 
Real-time Mode 



• Flush Pipeline on Halt 




13. Go to Debug and select the option connect by Pressing Reset Button on the board 



* /C671X AH_XDS510 Emuldtor/TMS320C671x_0 - C671x - Co 


File Edit View Project 


^ GEL Option Profile Tools DSP/BIOS 


i»t 


Breakpoints... 


r 




Assembly/Source Stepping 
Step Into 
Step Over 
Step Out 


k 




I 


Fll 
FIO 
Shift+Fll 


ffi I ^ 


IS3 








f mi 

+ C2 GEL files 


Run 
Halt 


F5 
Shift+F5 


W 


1 CJ Projects 


Animate 


Alt+F5 




Run Free 


Ctrl+F5 


{*> 




Low Power Run 


Ctrl+Alt+F5 


H* 




Run to Cursor 


Ctrl+FlO 






Set PC to Cursor 


Ctrl+Shift+FlO 


"{} 




Restart 


Ctrl+5hift+F5 




Go Main 


Ctrl+M 


{*} 




Multiple Operation... 




X 


Advanced Resets 


► 


& 




Reset CPU 


Ctrl+R 


^ 




Reset Emulator 
Halt on Reset 


Ctrl+5hift+R 






m 


| Connect Alt+C | 


m 


Restore Debug State 


m 


Thread Level Debugging 




•l 




Real-time Mode 






Rude Real-time Mode 




^ Flush Pipeline on Halt 
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14. Now Target is connected 



S-iD GEL files 
1 dl Projects 



The target is now connected 



S3. Disassembly 




0©( 


3 


00007C1C 


00000000 




A 1 


00007C20 


00000000 






00007C24 


00000090 






^>|D0007C28 


00008000 






00007C2C 


00000000 




fail 


00007C30 


00000000 






00007C34 


00000000 






00007C38 


00000000 






00007C3C 


00000000 






Bnnnn7r , 4n 


flFFFA41 n 




V 1 


< 




>_ 





ft j HALTED 



For Help, press Fl 
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15 Go to Debug and press Reset CPU that will initialize your memory. 





< /C671X AH_XDS510 EmuUitor/TMS320C671x_0 C671x Code Composer Studio [Disassembly] 


Q@K 


m File Edit View Project Q Q GEL Option Profile Tools DSP/BIOS Window Help 




_ 3 X 


^H 


Breakpoints.,. 

Assembly/Source Stepping 
5tep Into 
Step Over 
Step Out 


Fll 
F10 
Shift+Fll 


|fi **%*& 


mw 


t)= ■+)= | H)i «5i | O 




6 * 


li 




I d 


£l 


AT I Q i an 




?} 
?} 


f Files 
El- Q GEL files 
Cj Projects 


Run 

Halt 

Animate 

Run Free 

Low Power Run 

Run to Cursor 

Set PC to Cursor 

Restart 

Go Main 

Multiple Operation,.. 

Advanced Resets 


F5 

Shift+F5 

Alt+F5 

Ctrl+F5 

Ctrl+Alt+F5 

Ctrl+FlO 

Ctrl+5hift+F10 

Ctrl+Shift+F5 

Ctrl+M 


WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 
WOP 


* 
- 




Reset CPU Ctrl+R 


Reset Emulator 
Halt on Reset 


Ctrl+Shift+R 


m 
m 


Disconnect 
Restore Debug State 


Alt+C 


Thread Level Debugging 
Real-time Mode 
Rude Real-time Mode 




ml ,il 


s^J hUlZJ 


_>. 




^ o HALTED 


1/ Flush Pipeline on Halt 


Reset target CPU 







After Pressing Reset the CPU following window should appear in CCS. 



'4 /C671X AH_XDS510 Emulfltor/TMS320C671x_0 - C671x - Code Composer Studio - [Disassembly] 



File Edit View Project Debug GEL Option Profile Tools DSP/BIOS Window Help 



i » 



"TJI !tt%^ | # W 



til- 



~nr 



I # m m a | « ^ -4 



ov a igaElsHd 



JTiiel 

+ Q GEL files 
• CJ Projects 



-i^r 



o 
o 













< 



0000000 
0000004 
0000008 

ooooooc 

0000010 
0000014 
0000018 
000001C 
0000020 
0000024 
0000028 
000002C 
0000030 
0000034 
0000038 
000003C 
0000040 
0000044 
0000048 
000004C 
0000050 



FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 
FFFFFFFF 



.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 
.word 



Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 
Oxffff 



ffff a 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 

ffff 



v 



> 



■ r 



For Help, press Fl 



Now Board is ready for working real and non real time programs. 



60 



Experiment 1: Linear Convolution 
Procedure to create new Project: 

1. To create project, Go to Project and Select New. 



I /ChipMax_671 3/TMS320C6713J) C671x Code Composer Studio 



File Edit View 



dfi o 



ft 
{? 
{? 
ft 
{? 

^f 

la] 
HI 
IS] 



1 Files 

□ ge 



Debug GEL Option Profile Tools PBC DSP/BIOS V 



Open... 

Use External Makefile. . . 

Export to Makefile,,. 

Add Files to Project.., 

Save 

Hose 



I Jl* 


# ffl m A\ ■& 


%«fi|Q, 



Source Control 



Compile File 
Build 

Rebuild All 
Stop Build 
Build Clean 



Ctrl+F7 
F7 



Configurations.., 
Build Options.., 
File Specific Options... 
Function Level Options. . . 



Project Dependencies. . . 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 




2. Give project name and click on finish. 



Project Creation 



Project Name: \\ 

Location: | C:\CCStudio_v3.1 \MyPrqects\ _J 

Fiujeut Type. [Executable (.out) ^| 

Target |TMS320CG7XK "3J 



< Back 




Cancel 




Finsh 


Help 







(Note: Location must be c:\CCStudio_v3.3\MyProjects). 
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3. Click on File 



^>N 



ew L 



-^ Source File, To write the Source Code. 



Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Wir 




Open... 
Close 
Save 
Save As.. 
Save All 



Ctrl+O 
Ctrl+S 



DSP/BIOS Configuration 




Load Program... 
Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Ctrl+L 
Ctrl+Shift+L 



Load GEL... 
Data 

Workspace 
File I/O... 



Difference between files. 
Merge Files... 



Print... 



Recent Source Files 
Recent Workspaces 
Recent Program Files 
Recent Symbols 
Recent GEL Files 



Launch Setup 
Exit 

— n 



Ctrl+P 
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Aim : Linear Convolution of the two given sequences 

Mathematical Formula: 

The linear convolution of two continuous time signals x(t) and h(t) is defined by 



•i 



For discrete time signals x(n) and h(n), is defined by 



?? 






Where x(n) is the input signal and h(n) is the impulse response of the system. 

In linear convolution length of output sequence is, 
Length (y(n)) = length(x(n)) + length(h(n)) - 1 

Program: 

#include<stdio.h> 
main() 

{ 

int m=4; 
int n=4; 
int i=0,j; 
intx[10]={l,2,3,4,0,0,0,0}; 

inth[10]={l,2,3,4,0,0,0,0}; 

int *y; 

y=(int *)0x0000100; 

for(i=0;i<m+n-l;i++) 

{ 

y[i]=0; 

for(j=0;j<=i;j++) 

y[i]+= X [j]*h[i-j]; 

} 
for(i=0;i<m+n-l;i++) 

printf("%d ",y[i]) ; 

} 

Output: 

1, 4, 10, 20, 25, 24, 16. 



/*Lengh of i/p samples sequence*/ 
/*Lengh of impulse response Co-efficients */ 

/*Input Signal Samples*/ 
/*Impulse Response Co-efficients*/ 
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4. Enter the source code and save the file with ".C" extension. 



'"* /ChipMax_671 3/TMS320C671 3_0 - C671x - Code Composer Studio - Not Connected - [new inear.c] 



< v . File Edit View Project Debug SEL Option Profile Tools PBC DSP/BIOS Window Help 



iiGSHli^el^f 



~rj| tf & s & «%v sn? 



/*%■ 



| new linear, pjt _^J | D ebug 



?} 

?} 
H* 

"{} 
{+} 

m 
m 



5 Files 

+ CJ GEL files 

El-^ Projects 

B-jce£ new linear. pjt (Deb 

O Dependent Project 

O Documents 

Q D5P/BIOS Config 

| Generated Files 

CJ Include 

+ Q Libraries 

- ^ Source 

[i?l new inear.c 



"31 # b m A\ ■ ► -£ 



^ QiaiiBgd 



l#include<stdio .h> 

main( ) 

l{ int ni=4; ,--'*£ engnt of iSp samples sequence*,--' 
lint n=4; ,<-'*£engnt of impulse response Co -efficients *s 

int i=0,j ; 
lint x [ 10 ] ={ 1,2,3,4,0,0,0,0}; y'*Inpu t Signal Samples*,--' 
lint h[10]={ 1,2,3,4,0,0,0,0}; s'*Impulse Response Co-effi 
Y'-'*At tne end of input sequences pad *M* and W f no, of 
lint *y; 

y=(int *)0x0000100; 

for ( i =0 ; i <m+n - 1 ; i ++ ) 

I* 

y[i]=0; 

|for(j=0;j<=i;j++) 
ly[i]+=x[j]*h[i-j]; 

I 

for ( i = ; l <m+n - 1 ; i ++ ) 
printf ("%d\n",y[i] ) ; 
I 



cients*s' 
zero 's*,--' 



5. Right click on source, Select add files to project .. and Choose " .C " file Saved before. 



"4 /ChipMax_671 3/TMS320C671 3_0 - C671x - Code Composer Studio - Not Connected - [new inear.c] 



<^ File Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



i\ »Q I * Gtoel 



"TJ| H^ti IV 



\A%% 



| new linear.pjt 



^ Debug 



(? 
"ft 

ft 

2 

m 

m 
m 



Q iaiBBd 



f Files 
+ Q GEL files 
B ^ Projects 

- _^ new linear.pjt (Deb 

Pi Dependent Project 
Pi Documents 
£j DSP/BIOS Config 
r~l Generated Files 

I CJ Include 

1+ C3 Libraries 

+ _J 



«/ Allow Docking 
Hide 
Float In Main Window 



. M omo 



& 



M 3fe® BIQ 



l#include<stdio .h> 

maim ) 

l{ intm=4; ,>-'*£ engnt of i,--'p samples sequencers' 
lint n=4; ,<-'*Lenght of impulse response Co -efficients *,*-' 

int i=0, j ; 
lint x[10]={l,2,3,4,0,0,0,0}; s'*Input Signal Samples*,--' 
lint h[10]={l,2,3,4,0,0,0,0}; s'*Impulse Response Co-efficients*,--' 
Y'-'*At tne end of input sequences pad f M f and W no, of zero f s*,--' 
lint *y; 

y=(int *)0x0000100; 

for(i=0;i<m+n-l;i++) 



[i]=0; 

br(j=0;j<=i;j++) 

[i]+=x[j]*h[i-j]; 



Tor ( i = ; i <m+n - 1 ; i ++ ) 
printf ( "%d\n " ,y [i ] ) ; 
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6. Right Click on libraries and select add files to Project., and choose 
C:\CCStudio_v3.3\C6000\cgtools\lib\rts6700.1ib and click open. 



! KhipMax _671 3/TM5320C671 3_0 - C671 x - Code Composer Studio - N 



^File Edit View Project Debug GEL Option Profile Tools PBC D5P/BI05 U 



]£H|*%i 



"3 1 & 



new linear, pjt 



jj | Debug 



3#Ii 



SI* 1 QlHSId »? 



fTiies 
S-D GEL files 
g-Q Projects 

- ^ new linear.pjt (Deb 

_| Dependent Project 
■ jl Documents 
_| DSP/BIOS Config 

| Generated Files 

-HI Include 

ai 

+ [JSourcl 



#include<stdio.h> 

mainf) 

{ int m=4; /*Lenght 
lint n=4; /*Lenght of 
lint i=0,j; 

int x[10]={l,2,3,4,0 
lint h[10]={l,2,3,4,0 
u/*At the end of ispu 

int *y; 
|y=(int *)0x0000100; 
=0;i<m+n-l;i++) 



Add Files to Project 



yf Allow Docking 
Hide 
Float In Main Window 



=0; 
Kj=0;j<=i;j++) 
+=^[j]*h[i-j] 






? 

lfor(i=0;i<ra+n-l;i++) 
|printf("%d\n",y[i]]; 
} 



Look in: D lib 



~z\ <■ it] rrt & 



@]fascmath67x.lib Brts6400,lib 
g f astmath67xe , lib j§] rts6400_eh , lib 



grts6200,lib 
jjrts6200_eh.lib 
§]rts6200e.lib 
§]rts6200e_eh.lib 



grts6400e.lib 
§j rts6400e_e h,lib 

3 rts6700_eh,lib 



j§rts6700e,lib 
j§rts6700e_eh.lib 



File name: |rtsG700.lib Open I 

Files of type: | Object and Library Files (V; K .r) jj Cancel 



Help 



7. a)Go to Project to Compile . 

b) Go to Project to Build. 

c) Go to Project to Rebuild All. 



► File Edit View IBS Debug GEL Option Profile 



:-M 



lew linear.pjt 



S|tf 01 



+ CH GEL hi 

- Q Projec 
B- j nt 
_ 
_ 
_ 



New,,, 

Open,,, 

Use External Makefile,,, 

Export to Makefile,,, 

Add Files to Project,,, 

Save 



Source Control 




Euild :ir* 



Configurations,,, 
Build Options,,, 
File Specific Options,,, 
Function Level Options,,, 



Qfl+F7 
F7 



Project Dependencies,,, 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 



**i4m,xmmHVMUM*kumm**kTmMa 



^ File Edit View Project Debug GEL Option Profile 



il L* 



| new linear.pjt 



^|| Debug 



"311 < 



^ QiailBBd 



W 



§ Files 

H Q GEL files 

- ^ Projects 

- ^ new linear.pjt (Deb 

F~\ Dependent Project 
r~\ Documents 
Q D5P/BI05 Config 
CJ Generated Files 
] Q Include 

l+l-Hl Libraries 

B- ^ Source 



[new inear.c] "C:\CCStudio_v3.1\C60 



Compile Complete, 

Errors, Warnings, Remarks 
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8. Go to file and load program and load ".out" file into the board.. 



« /ChipMM_671 3/THS320C671 3_0 - C67 



Edit View Project Debug GEL Optior 



?| 

?: 
t 

2 



s 



New 
Open... 
Close 
Save 
Save As. 
Save All 



Ctrl+O 
Ctrl+S 



! 



Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Ctrl+Shift+L 



Load GEL... 
Data 

Workspace 
File I/O... 



Load Program 



Look in: Q Debug 
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Difference between files. 
Merge Files... 



Print. 
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Recent Program Files 
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Launch Setup 



Exit 

I I' W 



Ctrl+P 



File name: | new linear, out 



Open 



Files of type: ".out 



1 



Cancel 



Help 



A 
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9. Go to Debug and click on run to run the program. 




] File Edit View Project 
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10. Observe the output in output window. 
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1 1 . To see the Graph go to View and select time/frequency in the Graph, And give the 
correct Start address provided in the program, Display data can be taken as per user. 



1 /ChipMax_6713fTMS320C6713_0 - C671x - Code Composer Si 
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128 
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7 
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s 

On 
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OK 
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12. Green line is to choose the point, Value at the point can be seen (Highlighted by 
circle at the left corner). 



25.0 
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Experiment 2: Circular Convolution 
Procedure to create new Project: 

1. To create project, go to Project and Select New. 



!* /ChipMax_671 3/TMS320C671 3_0 - C671x - Code Composer Studio 



File Edit View | | Debug GEL Option Profile Tools PBC DSP/BIOS V- 
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■ 
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f Files 

s-CH eg 



Open... 

Use External Makefile. .. 

Export to Makefile... 

Add Files to Project... 

Save 

Close 



Source Control 



Compile File 
Build 

Rebuild All 
Stop Build 
Build Clean 



Ctrl+F7 
F7 



Configurations... 
Build Options... 
File Specific Options... 
Function Level Options. . . 



Project Dependencies. . . 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 



2. Give project name and click on finish. 



Project Creation 



Project Name: \\ 

Location: | C:\CCStudio_v3.1 \MyPrc|ects\ _J 

Fiujeut Type. [Executable (.out) t| 

Target |TMS320CG7XK "TJ 



< Back 



Finsh 



Cancel 



Help 



( Note: Location must be c:\CCStudio_v3.3\MyProjects ). 
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3. Click on File ' ^ New ' ^ Source File, to write the Source Code. 



Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Wir 




Open... 
Close 
Save 
Save As 
Save All 



Ctrl+S 



DSP/BIOS Configuration 

UT 



Load Program... 
Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Ctrl+L 
Ctrl+Shift+L 



Load GEL... 
Data 

Workspace 
File I/O... 



Difference between files. . . 
Merge Files... 



Print... 



Recent Source Files 
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Recent Program Files 
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Recent GEL Files 



Launch Setup 



Exit 



Ctrl+P 
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AIM: To implement circular convolution of two sequences. 

Circular Convolution: 

Let xi(n) and x 2 (n) are finite duration sequences both of length N with DFT's Xi(k) and 
X 2 (k). Convolution of two given sequences xi(n) and x 2 (n) is given by the equation, 

x 3 (n) = IDFT[X 3 (k)] 

X 3 (k) = Xi(k)X 2 (k) 

N-l 

x 3(n) = 2^ xi(m) x 2 ((n-m)) N 

m=0 



Program: 

#include<stdio.h> 

intm,n,x[30],h[30],y[30],ij,k,x2[30],a[30]; 
void main() 

{ 

printf("Enter the length of the first sequence:\n"); 

scanf("%d",&m); 

printf("Enter the first sequenceAn"); 

for(i=0;i<m;i++) 

scanf("%d",&x[i]); 

printf("Enter the length of the second sequenceAn"); 

scanf("%d",&n); 

printf("Enter the second sequenceAn"); 

for(j=0;j<n;j++) 

scanf("%d",&h|j]); 

/*If length of both sequences are not equal*/ 
if(m-n!=0) 

{ 

if(m>n) /* Pad the smaller sequence with zero*/ 

{ 

for(i=n;i<m;i++) 

h[i]=0; 
n=m; 

} 

for(i=m;i<n;i++) 

x[i]=0; 
m=n; 

} 

y[0]=0; 



73 



a[0]=h[0]; 

for(j=l;j<n;j++) /*folding h(n) to h(-n)*/ 

a[j]=h[n-j]; 

/*Circular convolution*/ 

for(i=0;i<n;i++) 
y[0]+=x[i]*a[i]; 
for(k=l;k<n;k++) 

{ 

y[k]=0; 
/* circular shift*/ 
for(j=l;j<n;j++) 
x2|j]=a[j-l]; 

x2[0]=a[n-l]; 
for(i=0;i<n;i++) 

{ 

a[i]=x2[i]; 

y[k]+=x[i]*x2[i]; 

} 
} 
/*displaying the result*/ 

printf(" the circular convolution is\n"); 
for(i=0;i<n;i++) 
printf("%d ",y[i]); 

} 



Output: 

enter the length of the first sequence 
4 

enter the length of the second sequence 
4 

enter the first sequence 
432 1 

enter the second sequence 
1111 

the circular convolution is 
10 10 10 10 
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4. Enter the source code and save the file with ".C" extension. 



(circular.pjt 



_^]| Debug 



"TJ|#^tt^ *%> 



?} 

{} 
«* 
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S-O GEL Files 
S -^ Projects 

- ^ circular.pjt (Debug] 

C3 Dependent Project 

I C3 Documents 

CJ DSP/BIOS Config 
C~\ Generated Files 
Pi Include 

Pi Source 



< 



l#include<stdio .h> 

lint m,n,x[30] ,h[30],y[30] ,i,j ,temp[30] ,k,x2 [30] ,a [30] ; 
void main( ) 

|{ 

int *y; 
ly=(int *)0x0000100; 

printf(" enter the length of the first sequence\n" ) ; 

scanf ( "Xd" ,&m) ; 
lprintf(" enter the length of the second sequence\n" ) ; 

scanf ( "Xd" ,&n) ; 
lprintf(" enter the first sequence\n" ) ; 
|for(i=0;i<m;i++) 

scanf ( "Xd" ,StX [i ] ) ; 

printf(" enter the second sequence\n" ) ; 
|for(j=0;j<n;j++) 

scanf ( "Xd" ,&h [ j ] ) ; 
lif (m-n ! =0) S*If length of hath sequences are not equal *<-' 

|{ 

lif (m>n) -'-■* Pad the smaller- sequence with z:er&*s' 

\{ 

I for ( i =n ; i <m ; i ++ ) 

h[i]=0; 
ln=m; 
I 
I for ( i =m ; i <n ; i ++ ) 

x[i]=0; 

m=n; 



5. Right click on source, Select add files to project .. and Choose " .C " file Saved before. 
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Add Files to Project 



Look in: | & circular 



0H 



"3 ■ S & M* 



(£) circular. C5_ 
£) Debug 



circ.c 



File name: 

Files of type: |c Source Files ( K .c; K .ccc) 



"3 



Open 



Cancel 



Help 



J 



6. Right Click on libraries and select add files to Project., and choose 
C:\CCStudio_v3.3\C6000\cgtools\lib\rts6700.1ib and click open. 



circular.pjt 



2J Debug 



iil 



AT QgAHQld 



?} 
}} 

1} 
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Files 

[]GEL files 
B Q Projects 

- ^ circular.pjt (Debug] 

[] Dependent Project 
!"■■•[] Documents 
I CJ D5P/BI05 Config 

| [] Generated Files 
| [] Include 



rr 



v' Allow Docking 
Hide 
Float In Main Window 






#include<stdio.h 
int m,n,x[30],h[ 
void iflin() 

{ 

int *y; 

y=(int *)0x00001 

printff 1 enter t 

scanf("&J",&m); 

printf(" enter t 

|scanf("Xd",&n); 

printf(" enter t 

for(i=0;i<i;i++ 

scanf(7d",&K[i 

printf(" enter t 

for(j=0;j<n;j++ 

scanf("Xd",&h[j 

if(m-n!=0) /*// 

{ 

if(m>n) /* Pad t, 

{ 

f rtT { i =Y\ ' 1 /ffl ' 1 4-4- \ 



Add Files to Project 



Look in: Q lib 



J + B & lh 



§]fast:math67x.lib @]rts6400.1ib 

§]fastmath67xe,lib §|rts6400_eh.lib 

§]rts6200.lib grts6400e.lib 

i|rts6200_eh,lib § rts6400e_e h.lib 

J|rts6200e,lib BIBS! 

||rts6200e_eh.lib @|rts6700 eh.lib 



i|rts6700e,lib 
§]rts6700e_eh,lib 



File name: 
Files of type: 



|rts6700.lib 



Type: LIB File 

Date Modified: 4/1/2005 2:24 

5ize: 1,84 MB 



PM 



Object and Library Files ( K .o x ; x .l x ) 



"3 



Open 



Cancel 



Help 
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7. a) Go to Project to Compile . 

b) Go to Project to Build. 

c) Go to Project to Rebuild All. 



^ File Edit View Project Debug GEL Option Profile Tools PBC D5P/BI05 Window Help 
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- £3 Projects 

- ^ circular.pjt (Debug] 

CJ Dependent Project 
CJ Documents 

I- □ D5P/BI05 Conhg 

I r~l Generated Files 

| CJ Include 

lt| -r~| Libraries 

+ CJ Source 



[circ .c] "C :\CCStudio_v3 . l\C6000\cgtools\bin\cl6x' 

Compile Complete, 

Errors, Warnings, Remarks. 



~~ | | | |\ Build /^Stdout/ 



^j • j HALTED: s/w breakpoint 



l#include<stdio .h> 

lint m,n,x[30],h[30],y[30],i,j , temp [30] ,k,x2[30],a[30]; 
void main i ) 

|{ 

lint *y; 

y=(int *)0x0000100; 

lprintf(" enter the length of the first sequenceXn" ) ; 
I scan f ( "%d" ,£m) ; 

lprintf(" enter the length of the second sequence\n" ) ; 
I scan f ( "%d" ,&n) ; 

lprintf(" enter the first sequence\n" ) ; 
|for(i=0;i<m;i++) 
I scan f ( "%d" ,&x [i ] ) ; 

lprintf(" enter the second sequenceXn" ) ; 
|for(j=0;j<n;j++) 
Iscanf ("%d",&h[j ] ) ; 

if(m-n!=0) s'#If longth of both s&quenc&s ctro not &qual*S 

|{ 

if (m>n) ,-'* Pad tie smaller sequence with zero*,-'' 

|{ 

| for ( i =n ; i <m ; i ++ ) 
lh[i]=0; 
n=m; 

I 

I for (i=m;i<n;i++) 
x[i]=0; / 

■ m=n ' / In case of 

any errors 
or 

warnings 
it will be 
displayed 
here 



-fr"C:/CCStudio v3 . 1/1 
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8. Go to file and load program and load ".out" file into the board.. 



Load Program 



Look in: Q Debug 



~3 - s &■ mm- 



File name: |circular.out 
Files of type: | K .out 



~Z 



Open 



Cancel 



Help 



J 



9. Go to Debug and click on run to run the program. 
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Reset CPU 
Restart 
Go Main 


Ctrl+R 

Ctrl+5hift+F5 

Ctrl+M 


Reset Emulator 
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Restore Debug State 


Ctrl+Shift+R 
Alt+C 


Enable Thread Level Debugging. . . 
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10. Enter the input data to calculate the circular convolution. 
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enter the length of the second sequence 
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in 
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[in circc 
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4 3 2 1 
enter the second sequence 



The corresponding output will be shown on the output window as shown below 
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enter the second sequence 
1111 

the circular convolution is 
10 10 10 10 
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f HALTED: s/w breakpoint 
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1 1 . To see the Graph go to View and select time/frequency in the Graph, and give the 
correct Start address provided in the program, Display data can be taken as per user. 
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Acquisition Buffer Size 


12B 






Index Increment 1 




Display Data Size 


4 






DSP Data Type 


32-bit signed integer 




Q-value 







Sampling Rate (Hz) 1 


Plot Data From 


Left to Right 




Left-shifted Data Display 


Yes 




Autoscale 


On 




DCValue 







Axes Display 


On 




Time Display Unit 


s 


Status Bar Display 


On 




Magnitude Display Scale 


Linear 






Data Plot Style 


Bar 




v 












OK Cancel | 


Help | 
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Experiment 3: N-Point DFT 
Procedure to create new Project: 

1. To create project, Go to Project and Select New. 



J /ChipMax_6713/TM5320C6713_0 - C671x - Code Composer Studio 



File Edit View I 1 Debug GEL Option Profile Tools PBC D5P/BIOS V 




B\<& 



?} 
0* 

?} 

0» 






H 



Files 

Qge 



Open... 

Use External Makefile. 

Export to Makefile... 

Add Files to Project... 

Save 

Close 



Source Control 



Compile File 
Build 

Rebuild All 
Stop Build 
Build Clean 



Ctrl+F7 
F7 



Configurations... 
Build Options... 
File Specific Options... 
Function Level Options. , 



Project Dependencies. . , 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 




2. Give project name and click on finish. 



Project Creation | ) 








Project Name: If 




Location: | C:\CCStudio_v3.1 \MyPrciects\ _J 




FiujeuL Type. | Executable (.out) ^| 




Target |TMS320CG7XK _;J 




Cancel Help 
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( Note: Location mus t be c:\CCStudio_v3.3\MyProjects ). 

3. Click on File ' ■> New ' ■> Source File, To write the Source Code. 



Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Wir 




Open... 
Close 
Save 
Save As.. 
Save All 



Ctrl+O 
Ctrl+S 



DSP/BIOS Configuration 




Load Program... 
Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Ctrl+L 
Ctrl+Shift+L 



Load GEL... 
Data 

Workspace 
File I/O... 



Difference between files. 
Merge Files... 



Print... 



Recent Source Files 
Recent Workspaces 
Recent Program Files 
Recent Symbols 
Recent GEL Files 



Launch Setup 
Exit 

— n 



Ctrl+P 
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AIM : TO COMPUTE N-POINT DFT OF A GIVEN SEQUENCE AND TO PLOT 
MAGNITUDE AND PHASE SPECTRUM. 

Discrete Fourier Transform: The Discrete Fourier Transform is a powerful 
computation tool which allows us to evaluate the Fourier Transform X(e )C0 ) on a digital 
computer or specially designed digital hardware. Since X(e JC0 ) is continuous and periodic, 
the DFT is obtained by sampling one period of the Fourier Transform at a finite number 
of frequency points. Apart from determining the frequency content of a signal, DFT is 
used to perform linear filtering operations in the frequency domain. 

The sequence of N complex numbers xo,..., x n-\ is transformed into the sequence of N 
complex numbers Xq, ..., X N -\ by the DFT according to the formula: 



N-l 

X(k) = Xx(n)e" j27rnk/N 
n=0 



k = 0,l, ....N-l 



Program: 

#include <stdio.h> 
#include <math.h> 
#defme N 4 



//number of data values 



float pi = 3.1416; 
short x[N] = {1,1,0,0}; 
float real[ 10], imaginary [10]; 
int index=0: 



//Output results 



void dft(short *x, short k, float *out) //DFT function 

{ 



float sumRe = 0; 
float sumlm = 0; 
int i = 0; 
float cs = 0; 
float sn = 0: 



//initialize real component 
//initialize imaginary component 

//initialize cosine component 
//initialize sine component 



for (i = 0; i < N; i++) 

{ 
cs = cos(2*pi*(k)*i/N); 

sn = sin(2*pi*(k)*i/N); 

sumRe = sumRe + x[i]*cs; 

sumlm = sumlm - x[i]*sn; 

} 

real [index] =sumRe ; 

imaginary[index]=sumlm; 

index++; 

printf("%f %f\n",sumRe,sumIm); 



//forN-pointDFT 



//real component 

//imaginary component 

//sum of real components 
//sum of imaginary components 

//sum of real components 
//sum of imaginary components 
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} 



void main() 

{ 

intj; 

for(j = 0;j<N;j++) 

{ 
dft(x,j,out); 

} 
} 
Output: 

2.000000 0.000000 
0.999996-1.000000 
0.000000 0.000007 
1.000011 1.000000 



//call DFT function 



4. Enter the source code and save the file with ".C" extension. 



■'< /C671 3 Device Cycle Accurate Simuldtor7TMS32QC671 3 - C671 3 (Simulator) - Code Composer Studio - [dftl .c *] 



^ File Edit View Project Debug GEL 

18 * t ^|f 



Option Profile Tools PEG D5P/BI05 Window Help 



3|tf ft 'fc*4 ft ?.£!#¥? til ->!§ | *!§ tfi | £1 H | , %% 



|DFT1.pjt 



^|| Debug 



^ # m EDDD P * 



?} 

{} 
{? 

"{} 
ft 

n 



^|^r Qi^[UHH£i|gi^ 



J Files 

+ CJ GEL files 

B-Q Projects 

E-£§ DFT 1. pjt (Debug) 

Cj Dependent Project 
Pi Documents 
CJ D5P/BI05 Config 
O Generated Files 
r~\ Include 
+ CJ Libraries 

+ r ~i HBB H 



I > 



-QjAT 



(float pi = 3. 1416; 
short x [N ] = {1,1,0,0}; x>7-^F£f7£? 
Ifloat out [2] = {0,0}; 



s'sinit Re and Im results 



Ivoid dft( short *x, short k, float *out) SSDFT function 

V 



float sumRe =0; 
float sumlm =0; 
int i = 0; 
float cs = 0; 
float sn = 0; 

for (l = 0; l < N; i++) 

{ 

cs = cos(2*pi*(k)*i/tSI) ; 
sn = sin(2*pi*(k)*i/H) ; 
sumRe = sumRe + x[i]*cs; 
sumlm = sumlm - x[i]*sn; 

} 
out[0] = sumRe; 
out [ 1] = sumlm; 
printf ("%£ %f\n" ,out [0] ,out[l] 1 



} 

Ivoid main f ) 

|{ 

I int j ; 



for (j = 0; j < N; j++) 

{ 

dft(x,j ,out) ; 

} 



} 



yyinit real component 
yyjnit imaginary.-" component 

Ssinit cosine component 
Ss'init sine component 

s's'for N-point DFT 

yy'real component 
SSimaginary componen t 
s's'sum of real components 
s's'sum of imaginary components 

s's'sum of real components 
s'ssum of imaginary components 



Sscall DFT function 
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5. Right click on source, Select add files to project .. and Choose " .C " file Saved before. 



<r/TMS320C6713-C6713(Sir 



► File Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



I GS q I A m 



|DFT1. F 



^ ||Pebug 



ll#iil • £ 



^ I i §3H]0[e] £ I £§^ 



la 



fife 

+ Q GEL files 
- £3 Projects 

B £§ DFTl.pjt (Debug) 

P~l Dependent Project 

| Documents 

CJ DSP/BIOS Config 
C~\ Generated Files 
Cj Include 

El d Libraries 

+ _l 



TJ| fi 1 ft * 1* ft *& | # HP %*\ 



[Float pi = 3.1416; 

I short x[N] = {1,1,0,0}; Ss'l -cycle 
Ifloat out [2] = {0,0}; 



/'/'init fie and 1m res: 



void dft( short *x, short k, float *out) 
|{ 

I float sumRe =0; 

I float sumlm =0; 

I int i = 0; 

I float cs = 0; 

I float sn = 0; 



i/ Allow Docking 
Hide 
Float In Main Window 



>r (i = 0; i < N; i++) 

{ 

cs = cos(2*pi* (k)*i^N) ; 
sn = sin(2*pi* (k)*i/W) ; 
sumRe = sumRe + x[i]*cs; 
sumlm = sumlm - x[i]*sn; 



'/'BFT function 

'/'init real component 
'/'init imaginary,? component 

in it cosine component 
in it sine component 

'for- N-point BFT 



} 

Ivo 



} 

out[0] = sumRe; 
out[l] = sumlm; 
printf {"%£ %f\n",out [0] ,out[l] ) ; 



/'/'real component 

/'/'imaginary componen t 
''--''sum of real components 
-'/'sum of imaginary components 

/'/'sum of real components 
/'/'sum of imaginary components 



void maim ) 

{ 
int j ; 

for (j = 0; j < N; 

{ 

dft(x,j ,out ) ; 



J++) 




'/'call BFT function 



Add files to the project 



"? 



f 



Ln 8, Col 1 



*>/ '£ 5 rr 
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6. Right Click on libraries and select add files to project and choose 
C:\CCStudio_v3.3\C6000\cgtools\lib\rts6700.1ib and click open. 



i Accurate Simulator/TMS320C671 3 - C671 3 (Simulator) Code Composer Studio [dftl .c *] 



► File Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



I GS q I A m 



|DFT1. F 



^ ||Pebug 



ll#iil • £ 



^ I i §3H]0[e] £ I £§^ 



la 



fife 

+ CJ GEL files 
- £3 Projects 

B ■■£§ DFTl.pjt (Debug) 

P~l Dependent Project 
| Documents 

CJ dsp/bios config 

C~\ Generated Files 
Cj Include 

l-m 

+ CJ Source 



tji fi 1 a * i* ft *.& | # up ^%| 



Ifloat pi = 3.1416; 

I short x[N] = {1,1,0,0}; Ss'l -cycle 
Ifloat out [2] = {0,0}; 



/'/'init fie and 1m resi 



<S Allow Docking 
Hide 
Float In Main Window 




void dft( short *x, short k, float *out) 
|{ 

I float sumRe =0; 

I float sumlm =0; 

I int i = 0; 

float cs = 0; 

iat sn = 0; 

(i = 0; i < N; i++) 
{ 

cs = cos(2*pi* (k)*i^N) ; 
sn = sin(2*pi* (k)*i/N) ; 
sumRe = sumRe + x[i]*cs; 
sumlm = sumlm - x[i]*sn; 



} 

Ivo 



out[0] = sumRe; 
out[l] = sumlm; 
printf {"%£ %f\n",out [0] ,out[l] ) ; 



void maim ) 

{ 
int j ; 

for (j = 0; j < N; 

{ 

dft(x,j ,out ) ; 



J++) 




Add files to the project 



'/'BFT function 

/'/'init real component 
/'/in it imaginary.' component 

/'/in it cosine component 
/'/in it sine component 

/'/'for N-point BFT 

/'/'real component 

/'/'imaginary componen t 
''/'sum of real components 
-'/'sum of imaginary components 

/'/'sum of real components 
/'/'sum of imaginary components 



'/'call DFT function 



"? 



f 



Ln 8, Col 1 



< ■;- e?- 



Add Files to Project 



Look in: Q lib 



"21 ^ m d° H 



!3]fasl:mal:h67x.lib 

3)fastmal:h67xe.lib 

||rl:s6200.lib 

i|rts6200_eh.lib 

9)rts6200e.lib 

3)rts6200e eh. lib 



||rr_s6400.lib 
||rr_s6400_eh.lib 
§|rt:s6400e.lib 
rts6400e eh. lib 



]rts6700e.lib 

]rts6700e eh. lib 



rts6700 eh .lib 






Type: LIB File 

Date Modified: 4/1/2005 2:24 PM 

Size: 1.34 MB 



Filename: |rts6700.lib 

Files of type: | Object and Library Files ( K .o K ; K .r 



Upen 



~3 



Cancel 



Help 



A 
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7. a) Go to Project to Compile . 

b) Go to Project to Build. 

c) Go to Project to Rebuild All. 




„*,,*, 



,,,.,,„ 



^ File Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS 



j£Q I 



DF11.P) 



ft 
(? 

ft 
(? 

H 

(1 
% 

9 
9 

H 

1 



I 



Files 
B-fiGELfil 
BQProjec 

i i 

i L 



+ _ 

- j 



i£Q| AMil*" 



New... 
Open.., 

Use External Makefile.. 
Export to Makefile... 
Add Files to Project.,. 
Save 
Close 



ft 



Source Control 




Stop 
Build Clean 



Configurations... 
Build Options.,. 
File Specific Options... 
Function Level Options., 



> 



Ctrl+F7 
F7 



Project Dependencies.,. 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 



DFTI.pjt 



jJI Debug 



3 #iasA * 



p|^ QISiBHii ii^&i 



f Files 

+ C3 GEL files 

B-Q Projects 

- ^j DFTI.pjt (Debug) 

Q Dependent Project 
Q Documents 
DSP/BIOS Config 
_| Generated Files 
Q Include 

+ C] Libraries 

H'Q Source 
] |I|dftl.c 



w 



y/BFT af N-pamt fr> 

l#include <stdio.h> 
l#include <math.h> 

void dft( short *x, 

#define N 4 

float pi = 3.1416; 
Ishort x[N] = {1,1,0 

float out[2] = {0,0 

jvoid dft( short *x, 

( 
float sumRe = 0; 

float sumlm = 0; 

int i =0; 

float cs = 0; 

float sn = 0; 

for (l = 0; l < N; 

{ 

cs = cos(2*pi 
sn = sin(2*pi 
sumRe = sumRe 

, rnniTtn = rntnTm 



out[ 

au± r 



It will shows 
warnings and 
errors if any 



[dftl.c] "C:\CCStudio_v3.1\C600 




s\bin\cl6s 



Compile Complete, 
Errors, Warnings, Remarks. 



| | | |\Md^Stdout/ 



8. Go to file and load program and load ".out" file into the board.. 



& /ChipMax_671 3/TMS320C671 3_0 - C67 



4jHga Edit View Project Debug GEL Optior 



New 
Open... 
Close 
Save 
Save As.. 
Save All 



Ctrl+O 
Ctrl+S 



Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Load GEL... 
Data 

Workspace 
File I/O... 



► 



Difference between files.. 
Merge Files... 



Print.. 



Recent Source Files 
Recent Workspaces 
Recent Program Files 
Recent Symbols 
Recent GEL Files 



Launch Setup 



Load Program 



Look in: & Debug 



§DFTl.out 



File name: 



Files of type: r.out 



~3 & B o- i^ 



Open 



~3 



Cancel 



Help 



9. Go to Debug and click on run to run the program. 




GEL Option Profile Ipols PBC DSP/BIOS Wi 



n i^y 



~3 



Breakpoints... 
Probe Points.,, 



a 



m 



new linear.pjt 



tt» 
{} 
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Step Into 
Step Over 
Step Out 



J Files 

+ CJ GEL files 

- £3 Projects 

^ new Mm 

CJ Depend 
icume 
'5P/BIC 
Generat 
Include 
F^ Libraries 
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Source 
\£\ new 




■3 

/Hall 

f Anin 



Halt 

Animate 

Run Free 

Run to Cursor 

Set PC to Cursor 

Multiple Operation... 



AssernblvjSource Stepping 



Reset CPU 
Restart 
Go Main 



Reset Emulator 
Disconnect 
Restore Debug State 



Fll 
F10 
Shift+Fll 



Shift+F5 

Alt+F5 

Ctrl+F5 

Ctrl+FlO 

Ctrl+Shift+FlO 



Ctrl+R 

Ctrl+Shift+F5 

Ctrl+M 



Ctrl+Shift+R 
Alt+C 



Enable Thread Level Debugging. . 

Real-time Mode 

Enable Rude Real-time Mode 



^ Flush Pipeline on Halt 

■UUULUAUH UUUUUb^ 
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89 



10. Observe the output in output window. 



ate Simulator/TMS320C671 3 - C671 3 (Si 



File Edit View Project Debug GEL Option Profile Tools PBC DbP/BIOS Window Help 
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^ IJDebug 
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+ CJ GEL files 
- ^ Projects 
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| Dependent Project 
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Q DSP/BIOS Config 

| Generated Files 
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H ■ CD Source 
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Output will 
shown here 



2.000000 0.000000 

0.999996 -1.000000 

0.000000 0.000007 

1.000011 1.000000 



~~ I I I \\ Build/, Stdoirt f~ 



^d *-> HALTED: s/w breakpoint 



For Help, press Fl 



1 1 . To see the Graph go to View and select time/frequency in the Graph, and give the 
correct Start address provided in the program, Display data can be taken as per user. 
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12. Green line is to choose the point, Value at the point can be seen (Highlighted by 
circle at the left corner). 
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(lj 0.999996) (1,-1) 



Lin Auto Scale 
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Experiment 4: FIR FILTER 

Procedure to create new Project: 

1. To create project, Go to Project and Select New. 



'4 /ChipMax_671 3/TMS320C671 3_0 C671x Code Composer Studio 



File Edit View | | Debug GEL Option Profile Tools PBC DSP/BIOS V 
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Add Files to Project... 

Save 

Close 



Source Control 



Compile File 
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Rebuild All 
Stop Build 
Build Clean 



Ctrl+F7 
F7 
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Build Options... 
File Specific Options.., 
Function Level Options. . . 



Project Dependencies. . . 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 




2. Give project name and click on finish. 



Project Creation 



Project Name: | 
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Fiujeut Type. [Executable (.out) 
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(Note: Location must be c:\CCStudio_v3.3\MyProjects). 

3. Click on File ' ■> New ' ■> Source File, to write the Source Code. 



Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Wir 




Open... 
Close 
Save 
Save As.. 
Save All 



Ctrl+O 
Ctrl+S 



DSP/BIOS Configuration 




Load Program... 
Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Ctrl+L 
Ctrl+Shift+L 



Load GEL... 
Data 

Workspace 
File I/O... 



Difference between files. 
Merge Files... 



Print... 



Recent Source Files 
Recent Workspaces 
Recent Program Files 
Recent Symbols 
Recent GEL Files 



Launch Setup 
Exit 

— n 



Ctrl+P 
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AIM : Realization of FIR filter (any type) to meet given specifications. The input can be a 
signal from Function Generator/Speech signal 

Finite Impulse Response (FIR) Filter: The FIR filters are of non-recursive type, 
whereby the present output sample is depending on the present input sample and previous 
input samples. 
The transfer function of a FIR causal filter is given by, 

N-l 

H(z) = X h(n)z n 

n=0 
Where h(n) is the impulse response of the filter. 
The Fourier transform of h(n) is 

N-l 

H(e )w ) = Z h (n)e _Jwn 

n=0 
In the design of FIR filters most commonly used approach is using windows. 
The desired frequency response H^e^) of a filter is periodic in frequency and can be 
expanded in Fourier series. The resultant series is given by, 71 

h d (n) = (1/2tt) J H(Or dw 

-71 
And known as Fourier coefficients having infinite length. One possible way of obtaining 
FIR filter is to truncate the infinite Fourier series at n = ± [(N-l)/2] 
Where N is the length of the desired sequence. 

The Fourier coefficients of the filter are modified by multiplying the infinite impulse 
response with a finite weighing sequence w(n) called a window. 

Where w(n) = w(-n) ± for |n| < [(N-l)/2] 
= for|n|>[(N-l)/2] 

After multiplying w(n) with hd(n), we get a finite duration sequence h(n) that satisfies 
the desired magnitude response, 

h(n) = h d (n) w(n) for |n| < [(N-l )/2] 
= for |n| > [(N-l)/2] 

The frequency response H(e JW ) of the filter can be obtained by convolution of H^e^) 
and W(e JW ) is given by, 
71 

HCO = (1/2ti) I H^e 1 ' 9 ) W^ 9 ) d9 

-71 
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Program: 

#include <stdio.h> 
#include "c6713dsk.h" 
#include "master.h" 
#include "aic23cfg.h" 
#include "dsk6713_aic23.h" 
#define FilLen 32 

// Low Pass Filter 

/*float FilCo[32]= { -0.0006,-0.0056,-0.0129,-0.0193,-0.0176,-0.0043, 0.0156, 0.0282, 
0.0190,-0.0129,-0.0483,-0.0545,-0.0065, 0.0927, 0.2064, 0.2822, 0.2822, 0.2064, 
0.0927,-0.0065,-0.0545,-0.0483,-0.0129, 0.0190, 0.0282, 0.0156,-0.0043,-0.0176, 
0.0193,-0.0129,-0.0056,-0.0006};*/ 

Note: //Generate the Coefficient using matlab filter design tool box // 
// High Pass Filter 

float FilCo[32] = { -0.0128,-0.0007,0.0068,0.0150, 0.0173, 0.0091,-0.0082,-0.0260,- 
0.0317,-0.0160,0.0200, 0.0611,0.0808,0.0481,-0.0803,-0.5875, 0.5875, 0.0803,-0.0481,- 
0.0808,-0.0611,-0.0200,0.0160,0.0317, 0.0260, 0.0082,-0.0091,-0.0173,-0.0150,-0.0068, 
0.0007,0.0128}; 

DSK6713_AIC23_Config config = {\ 

0x0017, /* DSK6713_AIC23_LEFTINVOL Left line input channel volume */ \ 
0x0017, /* 1 DSK6713_AIC23_RIGHTiNVOL Right line input channel volume *A 
OxOOff, /*2DSK6713_AIC23_LEFTHPVOL Left channel headphone volume */ \ 
OxOOff, /* 3 DSK6713_AIC23_RIGHTHPVOL Right channel headphone volume */ \ 
0x00 10, /* 4 DSK67 1 3AIC23ANAPATH Analog audio path control */ \ 
0x0000, /* 5 DSK6713AIC23DIGPATH Digital audio path control */ \ 
0x0000, /* 6 DSK6713_AIC23_POWERDOWN Power down control */ \ 

0x0043, /* 7 DSK6713AIC23DIGIF Digital audio interface format */ \ 
0x008c, /* 8 DSK67 1 3_AIC23_S AMPLERATE Sample rate control */ \ 

0x0001 /*9DSK6713_AIC23_DIGACT Digital interface activation */ \ 

}; 

DSK67 1 3_AIC23_CodecHandle hCodec; 

Int32 Line = 4; 

Int32 RInc = 3; 

Int32 LData; 

Int32 RData; 

Inf32 InChannel =0; 

Int32 OtChannel =0; 

Int32 LIndex = 0; 

Int32 RIndex = 0; 

Int32 LDataArr[128]; 

Int32 RDataArr[128]; 

Int32 LWP=0; 

Int32 LRP=100; 
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Int32 RWP=0; 
Int32 RRP=100; 
float LPrBuf[ 150]; 
float RPrBuff 150]; 
float FilOut= 0; 

main() { 

int i; 

LED=0x0; 

CSL_init(); 

// Filter Initialization 

for( i = ; i < 100; i++ ) LPrBufIi]=0.0; 

for( i = ; i < 100; i++ ) RPrBufIi]=0.0; 

// Initialize codec 

hCodec = DSK6713_AIC23_openCodec(0, &config); 

IRQ_globalEnable(); 

IRQ_enable(IRQ_EVT_RINT 1 ); 

IRQ_enable(IRQ_EVT_XINT 1 ); 

} 

void led() { 

static int cc = 1 ; 

LED = cc; 

// To Shift Pattern 

if(cc== 0x03) cc = 0x05; 

else if (cc == 0x05) cc = 0x06; 

else if (cc == 0x06) cc = 0x03; 

else cc = 0x03; 

} 

setpU200M(){ 

} 
voidread(){ 

int i = 0; 



if(InChannel==0){ 
InChannel =1; 
// Process Left Channel 

LData=MCBSP_read(DSK67 1 3AIC23DATAHANDLE); 
for( i = 0; i <= FilLen-2; i++) LPrBufli] = LPrBufli+l]; 
LPrBuf[FilLen-l] = (float) LData; 
FilOut = 0.0; 

for( i = ; i <= FilLen- 1 ; i++ ) FilOut += (FilCo[i] * LPrBufli]); 

LData = (int)(FilOut/16); 
LDataArr[LWP++] = LData; 
if(LWP>=128)LWP=0; 
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} 

else{ 



InChannel =0; 

// Process Right Channel 

RData=MCBSP_read(DSK67 1 3_AIC23_DATAHANDLE); 

RDataArr[RWP++] = RData; 

if(RWP>=128)RWP=0; 



} 



void write() { 

if(OtChannel==l){ 

MCBSP_write(DSK67 1 3_AIC23_DATAHANDLE,LDataArr[LRP++]); 
if(LRP>=128) LRP = 0; 

OtChannel = 0; 



} 
else{ 



MCBSP_write(DSK67 1 3_AIC23_DATAHANDLE,RDataArr[RRP++]); 
if(RRP>=128) RRP = 0; 

OtChannel = 1; 



4. Enter the source code and save the file with main.c extension. 



_671 3/TMS32QC671 3_0 - C671x 



<^ File Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



T 8 & i* ™ % u m w Usg -s= | «;i m I o n I a % % % | « 



|FIRG713M.pjt 



_r] | Debug 



-Tjl#aa^l ® %\^%. m% I 



m\&\Q lS3L10Hii 






5 Files 

+ □ GEL files 

El-Q Projects 

B ^ FIR6713M.pjt (Debi 

C~\ Dependent Project 

ft Documents 

_| D5P/BI05 Config 

| Generated Files 

■Si Include 
■ C~l Libraries 
El-O Source 
_£] boot.cmd 



"0®~ 



l#include <stdio.h> 

l#include "c6713dsk.h" 

l#include "master. h" 

l#include "aic23cfg.h" 

l#include "dskrj713_aic23 .h" 

l#include <std.h> 

l#include <swi .h> 

l#include <log.h> 

l#include <c6x.h> 

l#include <csl .h> 

l#include <csl_mcbsp .h> 

l*'* Length of sine wave table *-<-' 
l#defme SIHE_TABLE_SIZE 48 

\ Delta 

I/*floa t fil ter_ CoefffJ ={0. 00, 0. 00, 0. 00 ', 0, 00, O, 00 ', O. 00 ', 0. 00, O. 00, 
O, 00, 0. 00, 0, 00, O, 00, O. 00, 0. 00, 0. 00, 0. 00, 0. 00, O, 00, O, 00, 0. 00, Q. 00, 0, 00, 0. Oi 
0. 00, O. 00, O. 00, 0. 00, O. 00, 0. 00, 1, 00, 0. 00} ,■ *.■- 

\-'-'' Low Pass Fil tor 

\-*floa t fil ter_ CoefffJ ={2. 7152976341711 46e - 003, 2. 315319 (e - 004, - . 

\- 1.94171 ' 4 3226e- 002, 1,773008 Te-002, -. 

I - 6. 0594- 91e - 002, 8. 970, 2 L L 5 '£ J 48266s - 002, 3, 01 45 '29* 9i 4625e- 901, 4. 

1-6.0594- 91 e- 002, -6.7 71730e-002, 2. 11 3436751 136944e-903, , 
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5. Right click on source, Select add files to project and Choose main.c file Saved before. 



iMiumnmimi* 



ij...i.un,ij.un,ij.unij.i.ijj.i..i.i.yjiiiii.ii.i 

^ File Edit View Project Debug GEL Option Profile Tools DSP/BI05 Window Help 



l^H|llg 



"3 a & "U ™ % u, I m *? Usi *§ | «si m I n n | •* * ^ % | « 



I Test-Code-Del. pjt ~^[| Debug 



-TJ|#ffltt^|^ & • <& 



3 6V Q I13^E@d 



^Hbs 
+ Q GEL files 
- Q Projects 

- £§ Test-Code-Del.pjt ( 

| Dependent Project 

C] Documents 
CJ DSP/BIOS Config 
CJ Generated Files 
Q Include 

; CJ Libraries 

El CJ Source 



'DJ3" 



l#include <stdio.h> 
l#include "c6713dsk.h" 
l#include "master. h" 
l#include "aic23cfg.h" 
l#include "dsk6713_aic23 .h" 
Itfdefine FilLen 32 

l-';-'' Low Pass Filter 
Is* float Fi I Co [32 J 



yy High Pass Filter 
float FilCo[32] 



|DSK6713_AIC23_Config confi 
0x0017, ,<"* O DST'~ 



21: 


0x0017, 


/# 


/ DSK6713 


22: 


OxOOff, 


.--■* 


2 DSK6713 


23: 


OxOOff, 


y* 


3 DSK6713 


24: 


0x0010, 


,--'* 


4 DSK6713 


25: 


0x0000, 


,---* 


5 DSK6713 


26: 


0x0000, 


,-<-* 


6 DSK6713 



{ -O. 0006, -O. OO 9, -O, 0103, -O, 0176, -O, 0043. 

-O, 0483, -O. 0545, -O, 0065, O, 0927, 0, 2064, 0, 2322. 
-O. 0545, -O. 0483, -O. 0129, O, 0190, O. 0282, O, 0156. 
-O, 0056, -O. 0006}; *S 

{ -0.0128,-0.0007,0.0068,0.0150, 0.0173, 0.0091,-1 
0.0200, 0.0611,0.0808,0.0481,-0.0803,-0.5875, l 
-0.0611,-0.0200,0.0160,0.0317, 0.0260, 0.0082,-1 
0.0007, 0.0128}; 

3 = { \ 

AIC23_L£FTINVOL Left line input channel volume *. 
ALC23_RLGHTLmVL Right line input channel volume ■■ 
AIC23_L£FIHPVOL Left channel headphone volume */ 
ALC23_RLGH17fPVOL Right channel headphone volume *. 
AIC23_AM4PA17f Analog audio path control */ 
AIC23_fJlGPAlH Digital audio path control *.--' 




"] | | \\ Messages /^ Build 



DISCONNECTED (HALTED) 



7. Right Click on the source add other supporting .c files which configure the audio 
codec. 



/C671X AH XDS51QEn 



ir/TMS320C671x - C671x - Code Cc 



► 



i g* Q \i ^ m I * 



T est-Code-D el. pjt _^J | D ebug 



ll#iii|«« m-& 



tj| * a •<• ^» <% Vb | « 'us 1 ( t^4|i 



& 6V Q lS3[|]0id 



T 


Files 


E- 


£i GEL files 


- Q Projects 




- J^ Test-Code-Del.pjt ( 




:■■ CJ Dependent Project 




f~~l Documents 




■ Q DSP/BIOS Config 




1 Generated Files 




:■■ CD Include 




:■■ -CD Libraries 




B £3 Source 








J3 aic23.c 








JU aic23cfg_c.c 








3H aic23cfg_csl.c 






J) dsk6713 aic23 








1) dsk6713 aic2J 








1) dsk6713 aic23 








^ main.c 






1 pll.c 


< 


Nil > 



-Djdf 



1: ^include <stdio.h> 



Look in: Test-Coda Del 



~3 <- E EJ Hh 



QTest-Code-Del.C5_ 

13 aic23 

[3 aic23cfg_c 

aic23cfg_csl 

dsk6?13_aic23_opencodec 
(3 dsk6713_aic23_registers 



(3 dsk6713_aic23_write 

[3 main 

©Pll 



r 



Files of type: |c Source Files ("-c;".ccc;".c1 6; K . c32) 



"31 



Help 



UxUUl/, 
OxOOff, 
OxOOff, 
0x0010, 
0x0000, 
0x0000, 



1 DSKh/L hUuHTlNVlh 

2 BSK6713_AIC: 

3 DSK6713_AIC 

4 I)SK6713_AIC 

5 BSK6713_AIC 

6 £JSK6713_AIC 



, -O. 0129, -O. 0193, -O. 0176, -O. 0043. 
0,0065, 0,0927, 0.2064, 0.2322. 
0.0129, 0.0190, 0.0282, 0.0156. 



,0.0068,0.0150, 0.0173, 0.0091,-1 
,0.0808,0.0481,-0.0803,-0.5875, l 
,0.0160,0.0317, 0.0260, 0.0082,-1 

}; 



r 1 



Left line input channel volume 



flight line input channel volume ■■ 
3_LEFTHPVOL Left channel headphone volume */ 
^23_RIGBTHPVOL Right channel headphone volume *- 
WAIN Analog audio path control */ 
3_I)IGPAJN Digital audio path control »S 

N Power down control *S J 
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8. Right Click on the source add other supporting .asm files which configure the audio 
codec, i.e. boot.asm 



'/C671X AH XDS510En 



ir™S32QC671x 0-C671x-C 



A File Edit View Project Debug 



E 

DSP/BIOS Window help 



m& y & %J 



JTest-Code-Del.pjt ^|| Debug 



H#iiiH)« •■£ 



?} 

?> 
I? 

m 
n 
ft 

m 





T||iff ai*V%?-i|#iff « >^%|« 



Pov aiais@fi 



f Files 


+. U GEL files 


B €3 Projects 


B £z£ Test-Code-Del.pjt ( 




| Dependent Project 




• Q Documents 




• Q DSP/BIOS Config 




■ CJ Generated Files 




■ O Include 




■-Q3 


Libraries 


- lA 


Source 






|j aic23,c 






J) aic23cfg_c.c 






H) aic23cfg_csl.c 






jK] boot.asm 






^ dsk6713 aic2: 






D dsk6713 aic2c 






^ dskj6713 aic2S 






ij main.c 




® pn.c 


< 



1: ^include <stdio.h> 



Look in: Test-Lod-;: D el 



~z\ <- is a* n~ 



(^ Test-Code-Del, C5_ 



fbooT 



Open || 



Files of type: JAsm Source Files ( K .a K ) 



t Cancel 



UxUUl/, 
OxOOff , 
OxOOff , 
0x0010, 
0x0000, 
0x0000, 



1 lfcKb/13_AlL 

2 r/SK6713_AIC 

3 DSK6713_AIC 

4 BSK6713_AIC 

5 BSK6713_AIC 

6 BSK6713_AIC 



" 



, -Q. 0129, -Q. 0193, -O. 0176, -Q. 0043. 
',-0.0065, 0.0927, 0.2064, 0.2022. 
',-0.0129, 0.0190, 0.0282, 0.0136. 



,0.0068,0.0150, 0.0173, 0.0091, 
,0.0808,0.0481,-0.0803,-0.5875, 
,0.0160,0.0317, 0.0260, 0.0082, 

}; 



■A. 



eft line input, channel volume * 



'2J_A'lttfll.: ■■■■ 

3_LEF7HP¥QL Left channel headphone volume #s' 

?23_J?flWFHP¥Q£ Right channel headphone volume *, 
3_MI%FA1H Analog audio path control */ 
3_BIGFA1H Digital audio path control »/ 
3_PQWEfiDOm Power down control *,--' ^J 




9. Right Click on the source add other supporting .asm files which configure the audio 
codec, i.e. boot.asm 



; 



File Edit View Project Debug GEL Option 



Pi ol He >/■.' Melo 



~t\\ & ?** -"ft *.& m *? | t ^ 4 1 •' 



E » S I I 



I Test-Code-Del. pit ~^"1| Debug 



"311 $ ffl m 



ov aniBgd 



ft Files 
S LJ GEL Files 
B £3 Projects 

- _^ Test-Code-Del.pjt ( 

I "1 Depend* nt i o 

| Documents 

_| D5P/BIOS Config 

I Generate ,i_ - 

O Include 

□ Libraries 

;-: £3 Source 

Q aic23.c 

J aic23c .,_..«. 

_if] boot.asm 

JO dsk£713_aic2 

j*T dskj6713_aic2: 

_; df t s-ij_=n:;: 

M main.c " 
Hpll.c 
i:t1 boot.crnd 



^ 



#include <stdio.h> 




Files of type: IQ 



■3 



Help 



, -tf. 9129, -O. 11193, -8. IU76, -9. 9943. 
O.LW65, 0. 11927, 8.2864, 11.2822. 
9. 9129, 9. 91911, 11.9282, 9.9188. 



,0.0063,0.0150, 


0.0173, 


0.0091,- 


,0.0803,0.0481, 


-0.0803, 


-0.5875, 


,0.0160,0.0317, 

}; 


0.0260, 


0.0082,- 



ef/'t 11ns put ahsnnsl vol urns *. 
'j?i£rnt line; input channel volume ■■ 
Loft ana. olume *<-'* 



LIxLiUl/ 1 , .-■'* 1 l&'A : 67llJ_Jj>"l "."<"'" ' 

OxOOff, ,-"* 2 i> r 7HFVOL 

OxOOff, .'<■* 3 'At cnannel headphi 

k 00 10, ---'* -;■ \ 'a si , trol *.-■' 

Gi;0u00, ."'* 5 Blffltal audio patn oontjrol *- 

, 



^ 



"1 



| | | |\ Messages /.Build/ 



HLiJ_ 
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1 0. Right Click DSP/BIOS Config on the add aic23 .tcf file which configure the 
DSP/BIOS. 



'/C671X AH XDS510Err 



TTHS320C671x Q-C671x 



i^ 


Help 








iu?u - 1| jia 'W^ 


mw 


»*|f 




| Test-Code-Del. pit _^]| Debug jj | <§* M -M S\ 


• £ 









r^|6V I o issiDBii I 



f Files 

+ _| GEL files 

- _j Projects 

B ■■■£ Test-Code-Del.pjt ( 

•■ O Dependent Project 
■ Cj Documents 
:-: £3 DSP/BIOS Config 

r aic23,tcf 
* Q Generated Files 
I--0 Include 
_| Libraries 
B Q Source 

1] aic23.c 
_£] aic23cfg_csl.c 
£\ boot. asm 
_£] dsk6713_aic2; 
£] dsk6713_aic2; 
£] dsk6713_aic23 
1 main.c 
H pll.c 
L£] boot.cmd 



^include <stdio.h> 



Look in: | £3 Test-Code-Del 



~3 <- IS ep n~ 



(fl| Test-Code-Del. CS_ 
: faic23 



I 



File i,i. re Configuration File f/\ tcf) 



"3 



UxUUl/, 
DxQOff, 
DxOOff, 
0x0010, 
0x0000, 
0x0000, 



'hyi3_AJL 

■ 2 DSK6713_AIC 

■ 3 DSK6713_AIC 

■ 4 DSK6713_AIC 

■ 5 BSK6713_AIC 

■ 6 BSK6713_AIC 



*k 



' 



-O. 0129, -0. 17193, -O. 0176, ■ 
-0,0065, 0.0927, 0,2064, 
-0.0129, 0.0190, 0.0282, 



■0. 0043. 
0, 2822. 
0, 0156. 



,0.0068,0.0150, 0.0173, 0.0091,- 
,0.0808,0.0481,-0.0803,-0.5875, 
,0.0160,0.0317, 0.0260, 0.0082,- 



'0ft line ., pel volume *, 

'23_A'lLiHl . '•• ' "Fight line input channel volume 
Left channel headphone volume *,■■' 



?_LEF7HPVOL 

?_FlGf£LHFVOL Fight ohannel headphone volume 

g audio path control *,--' 
3_BIGFA1H Digital audio path control *, 



J 




11. 



Right Click on Libraries add csl6713.1ib from 

C:CCStudio v3.3/C6000/csl/lib/csl6713 file which include the processor libraries. 



fy File Edit View Project Debug GEL ( 



ii & 



"3 1 ^.%^i#v? u %**\* 



> 



t || Debug 



"TJ|#ffiif» B%\ 



3 j 6V Qi®3L10®d 



f Files ' 

:E CH GEL files 
B'Q Projects 

- ^ Test-Code-Del.p 

• C~J Dependent Pro 

• Cj Documents 

- j DSP/BIOS Conf 
L-gP aic23.tcf 

- _J Genet al . il 

; [I) aic23cfg.se | 

1 LI aic23cfg_c 

| ■■ Q Include 

E! ^ Libraries 

■ TJ csl6713.lib 

S Si Source 

| ||j) aic23.c 

j Qj) aic23cfg_c: J, 

! H boot.asm 

| [D dsk6713_a 

; LI dsk6713_a. 

i~[H dsk6713_a 

# pll.c 
: +n boot.crmd 



"DJ3" 



1: #include <stdia.h> 



Look in: j Q 

Hcsl6413 
Hcsl6413e 
Hcsl6414 
Hcsl6414e 
Hcsl6415 
Hcsl6415e 



~3 «- & c* ^ 



|fflcsl6416 
|fflcsl6416e 
Hcsl6418 
Hcsl6418e 
Hcsl6701 
|fflcsl6701e 



gcs!6711 

gcsl6711C 

gcsl6711Ce 

Hcsl6711e 

Ocsl6712 

Hcsl6712C 



SJcsl6712Ce 
fflcsl6712e 



Jijcs!67i3e 
gcslDA610 
HcslDA610e 



File name: |cslG713 

Files of type: Object and Library Files [".o";M") 



~H 



Help 



UxUUlV, 

QxQOff, 
OxOOff, 
0x0010, 
OirOOOO, 
GkOOOU, 



, -O. 0129, -O. 0193, -O. 0176, -O. 0043. 

,-0,0065, 0.0927, 0.2064, 0.2822. 

0,0129, 0.0190, 0,0282, 0.0156. 



.0066, 
.0808, 
.0160, 



.0150, 
.0481, 
.0317, 



.0173, 
.0803, 

.0260, 



.0091,- 

.5875, 

.0032,- 



- j mbyyj_As:.-:'- — - — ■ — ' -i wjum* 

dphone volume *---' 
; 3 1/SX671. Sight channel headphones volume *, 

■ i _■"„"/". AFA1H Analog audio path control *=.<" 

: 5 DSK671 GPA7H Digital audio path control *.--' 

■ 6 BSKtl - -. o rol */ 



'■ 



■J 
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12 Go to Project- -> Build option -> Preprocessor^ set CHIP6713 as shown 
below. 



' /C671X AH_XDS510 Emulalor/TMS32GC671x_G - C671x - Code Composer Studio - Not Connected - [csL_chiphal.h] 



^ File Edit View Project Debug GEL Option Profile Tools DSP/BIOS Window 


Help 








11 fi? g | A -% B * || J| g ft^V^^|#V? 


» -Si | «Sl &i 1 S H 1 * r * £ % 1 < 






»-$ 








|Test-Code-Del.pit *j| Debug _*\ | & M M eI | « * 





ov q issmEStfi 



f? 

2 



§ Files 

+, CJ GEL hies 

R^ Projects 

B -£ Test-Code-Del.p 

! £] Dependent Pro 

I £] Documents 

B Si dsp/bios conf 

i--gpi aic23.tcf 
B £3 Generated File: 
d) aic23cfg.se 
[^1 aic23cfg_c 
r~l Include 
E! Q Libraries 

H csl6713.lib 
- £3 Source 

I] aic23.c 
I] aic23cfg_c: 
: £] boot, asm 
I] dsk6713_a 
Jj] dsk6713_a 
Jj] dsk6713_a 
|j main.c 
U pll.c 
liKl boot.cmd 

> 



25 
251 

25; 

251 
25' 
261 
26 
26; 
26: 
26- 
26! 
261 
26; 
261 
26- 
271 
27 

27; 
27; 

27- 
27! 
271 

27; 

271 
27' 
231 



; for Test -Code -De l.pjt (Debug) 



General Compiler | Linker j DspBiosBuilder | Link Order ] 



- r ■ 5-fr":| -i >' Debug" -d"_DEBUG" -mv6700 



Category: 







Diagnostics 



0: 



Preprocessor 
Include Search Path (- 
Pre-Define Symbol [-dj: |CHIP_6713J 
Undefine Symbol (-u): | 
• i\ cessing: JNone 

I - Continue with Compilation (-ppa) 



~3 




- 1 



sre XXXX is chip number, i.e. 6; 



c3,c9,cl0,cll,cl2,cl3,cl4,cl5, 



1 fatal error detected in the compilation of 
Compilation terminated. 

>> Compilation failure 



I j 1 |\ Messages \ Build f 



-> 
J 1 



7. 7. a) Go to Project to Compile. 

b) Go to Project to Build. 

c) Go to Project to Rebuild All. 
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* /Chipmax_671 3/THS320C671 3_0 - C671x - Code Composer Studio - [main.c] 



fy File Edit View 



Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



£ sSHU ' 



|FIR6713M.pit 



£3 &? a I 






1> 

m 
m 



New... 

Open... 

Use External Makefile. . . 

Export to Makefile... 

Add Files to Project,,. 

Save 

Close 




Ctrl+F7 
F7 



* 



Configurations... 
Build Options,,. 
File Specific Options,,. 
Function Level Options.. 



Project Dependencies, , . 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 



D File View j/^ Bookmarks j 



"TJ| 8 & fl k ™%?4 #*? 



m m a I * 3U^9 



> • B | Q 



* 



;5i -Si 



[^include <stdio.h> ~ 

#include "c6713dsk .h" 

#include "master. h" 

#include "aic23cfg.h" 

#include "dsk6713_aic23 .h" 

#include <std.h> 

^include <swi .h> 

^include <log.h> 

#include <c6x.h> 

#include <csl .h> 

#include <csl_mcbsp .h> 

,'''* Length of sine wave table *s 
^define SINE_TABLE_SIZE 48 

yy' Delta 

y*floa t fil ter_ Coefffj ={0, 00 ', tf . tftf, tf , < 
O, 00 ', A Z.7/.7, tf , /_7/_7, tf, tftf, /_7, Z7/.7, tf , T.7/.7, tf, /_7/_7, Q. t 
Q. 00 ', tf. tftf, tf, tftf, tf, tftf, tf, tftf, /_7, tftf, 1, 00 \ tf. * 

^-'' Zcw Pass Filter 

l^ float filter_Coeff[J ={2, 7152976341, 
1-1, 941705440790678e-002 r -1, 72909821 8t 
\-6, 059440700570791e - 002 ', tf , 9703132564- 
1-6, 059440700570791e-002 r -6, 788468549, 
1-1, 941705440790678e-002 r 2, 7341341665 

I V? -7 / C 2 <7 -7/." ^ / "7 7 / .^/.,V> _ O.O. 1 ? ■ * .-- 



«] 



D File View ) /4 Bookmarks 



I''.-' low Jr'ass ruzer 

/* float filtor_Cooff[] ={2, 71 52976341 71 146o- 883^ 2, 31581980294292 
- 1 . 941786448798678& - 882 r - 1 . 72989821 884 3226o - 882 ', 1 . 7738887385686/ 



I - 6. 85 i ?lo - 882 r 8. 97831325 ' 

I - 6. 859448788578791& - 882 r - 6. 7884685 
\-l. 941 786 4 487986 78& - 882 r 2, 73413416 



[Linking. . . ] "C:\CCStudio_v3. l\C6000\cgtools\bin\cl6x^- 
<Linking> 

Build Complete, 

Errors, Warnings, Remarks. 



It will shows 
warnings and 
errors if any 



—y374625 
T113694 
U54354 



j | | |\ Messages }\ Build f~ 

% m|HALTED 
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8. Go to file and load program and load ".out" file into the board. 



3 /Chipmax_6713/TMS320C6713_0 - C671x - Code Composer Studio - [n 



Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



i 



rt 

r 
K 
C 



New 
Open,,, 
Close 
Save 
Save As.. 
Save All 



Ctrl+O 
Ctrl+S 



| S| iMT fe ® [ Ej I Q> 



Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Ctrl+Shift+L 



Load GEL... 
Data 

Workspace 
File I/O... 



Difference between files, . , 
Merge Files,., 



Print... 



Ctrl+P 



Recent Source Files 
Recent Workspaces 
Recent Program Files 
Recent Symbols 
Recent GEL Files 



Launch Setup 



Exit 



l< 



▼ I | h M M Wh T4 ■■+ ^ nf- 



~T]| # ffl £3 ol « l|M* ► £ 



> 



r(= .(= (z r+Cz 



HI/*' 



^include 
^include 
#include 
^include 
^include 
^include 
^include 
#include 
^include 
^include 
^include 



<stdio .h> 

"c6713dsk.h" 

"master . h" 

"aic23cfg .h" 

ll dsk6713_aic23.h l 

<std.h> 

<swi .h> 

<log.h> 

<c6x.h> 

<csl .h> 

<csl_mcbsp .h> 



s* Length of sine wave table *>-'' 
#define SINE_TABLE_SIZE 48 

V Delta 

Y''*floa t 111 ter_ Coefff] ={0, 00 ', O. 00 } * O. 00 } * tf , tftf, O. 00 '* tf , /_7/J 
|/_7, tftf, /_7, tfft. /_7, tftf, /_7, tftf, tf, tftf, tf , tftf, tf , tftf, /_7, 00* O. 00* O. 00 } * tf, /_7/J 
Itf, tftf, /_7, tftf, tf. tftf, 7.7, 00* O. OO r O. 00* 1. 00* O. 00}; *s' 

\'''>-'' Low Pass Filter- 

Y"* float filter_Coeff[J ={2, 715297634171 146e-003* 2,315 
\-l, 941706440790678e-002* -1, 729098218843226e-002* 1, 775 
\-6. 059440700570791e-002* 8, 97031 3256 4 48266e- 002* 3, 0145 
\-6. 059440700570791e-002* -6, 788468549771 730e-002* 2, 115 
1-1. 941706440790678e-002* 2, 7341 341 66535232e- 003* 1, 2075 



Load Program 



Look in: | fcJH> Debug 






E C*" Unn- 



File name: |aic23.out 

Files of type: | K .out 



Open 



Cancel 



Help 



^w«* 



Ud 
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9. Go to Debug and click on run to run the program. 




] File Edit View Project 



i\ & 



new linear, pjt 



"3 






aiasi 



§ Files 
El-O GEL files 
B--4E3 Projects 
- ^ new line 




GEL Option Profile Tools PBC DSP/BIOS Wi 



Breakpoints... 
Probe Points... 



Step Into 
Step Over 
Step Out 



Halt 

Animate 

Run Free 

Run to Cursor 

Set PC to Cursor 

Mu[tiple Operation 



AssemblyjSource Stepping 



Reset CPU 
Restart 
Go Main 



Reset Emulator 
Disconnect 
Restore Debug State 



Fll 
F10 
Shift+Fll 




5hift+F5 

Alt+F5 

Ctrl+F5 

Ctrl+FlO 

Ctrl+5hift+F10 



Ctrl+R 

Ctrl+Shift+FS 

Ctrl+M 



Ctrl+Shift+R 
Alt+C 



Enable Thread Level Debugging. 

Real-time Mode 

Enable Rude Real-time Mode 



^ Flush Pipeline on Halt 

I UUUUMUH UUUUUb^ 



ni 



10. Output can be seen on CRO with the filtering effect. 
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Experiment 5: To plot spectrum of cosine signal plus a signal from 

function generator 
Procedure to create new Project: 



1. To create project, Go to Project and Select New. 



'4 /ChipMax_6713/TMS32QC6713_0 - C671x - Code Composer Studio 



File Edit View 



i\ & 



r 



SI 6V 



?} 

?} 
-ft 

m 



f Files 



Debug GEL Option Profile Tools PBC D5P/BI05 V 



Open... 

Use External Makefile. . 

Export to Makefile,,, 

Add Files to Project... 

Save 



Source Control 



Compile File 
Build 

Rebuild All 
Stop Build 
Build Clean 



Ctrl+F7 
F7 



Configurations... 
Build Options.., 
File Specific Options... 
Function Level Options, , , 



Project Dependencies, . , 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 




2. Give project name and click on finish. 



Project Creation 









Project Name: [| 








Location: | C: \CCS tudio_v3. 1 \MyPrcjects\ 


J 






Pi ujeul Type. | E xecutable (. out) 


"31 






Target TMS320C67>O< 


"31 






< Back Finsh 


Cancel 


Help 





(Note: Location must be c:\CCStudio_v3.3\MyProjects). 
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3. Click on File ' ■> New ' ■> Source File, to write the Source Code. 



J Edit 


View 


Project 


Debug GEL 


Option 


Profile Tools 


PBC 


DSP/BIOS Wir 


hmihh ^■^1 0V3FV WBEm^M 


Open... 






Ctrl+O 
Ctrl+S 




D5P/BI05 Configuration. . . 


Close 


J|fiW 


|±*| 


IS eB=i 1 









Save As.. 
Save All 



Load Program... 
Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Load GEL.. 



Workspace 
File I/O... 



Difference between files.. 
Merge Files... 



Print.. 



Recent Source Files 
Recent Workspaces 
Recent Program Files 
Recent Symbols 
Recent GEL Files 



Launch Setup 




AIM: To add a cosine signal to a signal from a function generator/Microphone, and 
to plot the spectrum. 

#include <stdio.h> 
#include <csl.h> 
#include"c6713dsk.h" 
#include "master.h" 
#include "aic23cfg.h" 
#include "dsk6713_aic23.h" 

/* Length of sine wave table */ 
#defme SINETABLESIZE 48 

// Low Pass Filter 

float filter_Coeff[] ={2.715297634171 146e-003, 2.315819802942924e-004, 

1.244493581373643e-002, 7.244364917221401e-003J.207341716154354e-002, 

2.734134166585232e-003,-1.941706440790678e-002,-1.729098218843226e- 

002J.773008730568675e-002,4.091495174059349e-002,2.1 13436751 136944e-003, - 

6.788468549771730e-002 9 -6.059440700570791e-002, 8.970313256448266e- 
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002,3.014572949374625e-001,4.019009454299968e-001,3.014572949374625e-001, 
8.970313256448266e-002,-6.059440700570791e-002,-6.788468549771730e- 
002,2.1 13436751 136944e-003, 4.091495174059349e-002,1.773008730568675e-002, - 
1 .7290982 1 8843226e-002,- 1 .94 1 706440790678e-002, 
2.734134166585232e003,1.207341716154354e-002,7.244364917221401e- 
003,1.244493581373643e-002,2.315819802942924e-004,2.715297634171146e-003}; 

// High Pass Filter 

/*floatfilter_Coeff[]={3.294316420702696e-004,3.800020076486443e- 
003,9.822200806739014e-003,1.517265313889167e-002, 1.323547007544908e- 
002,2.635896986048919e-004,-1.808215737734512e-002,-2.666833013269758e-002,- 
1.155354962270025e-002,2.448211866656400e-002,5.534101055783895e 
002,4.424359087198896e-002,-2.922329551555757e-002,-1.473332022689261e-001,- 
2.574625659073934e-00 1 ,6.976203 1 09 125493e-00 1 , -2.574625659073934e-00 1 , 
1 .47333202268926 1 e-00 1 ,-2.92232955 1 555757e-002,4.424359087 1 98896e-002, 
5.534101055783895e-002,2.448211866656400e-002,-1.155354962270025e-002, 
2.666833013269758e-002, -1.808215737734512e-002,2.635896986048919e 
004,1.323547007544908e-002,1.517265313889167e-002,9.822200806739014e- 
003,3.800020076486443e-003,3.294316420702696e-004};*/ 

// Pre-generated Co-sine wave data, 16-bit signed samples 

Int32 sinetable[SINE_TABLE_SIZE] = { 0x00000000, 0x000010b4, 0x00002120, 
0x000030fb, 0x00003fff, 0x00004dea, 0x00005a81, 0x0000658b,0x00006ed8, 
0x0000763f, 0x00007bal, 0x00007ee5, 0x00007ffd, 0x00007ee5, 0x00007bal, 
0x000076ef, 0x00006ed8, 0x0000658b, 0x00005a81, 0x00004dea, 0x00003fff, 
0x000030fb, 0x00002120, 0x0000 10b4, 0x00000000, 0xffffef4c, OxffffdeeO, 0xffffcf06, 
0xffffc002, 0xffffb216, 0xffffa57f, 0xffff9a75, 0xffff9128, 0xffff89cl, 0xffff845f, 
OxffffSllb, 0xffff8002, OxffffSllb, 0xffff845f, 0xffff89cl,0xffff9128, 0xffff9a76, 
0xffffa57f, 0xffffb216, 0xffffc002, 0xffffcf06, OxffffdeeO, 0xffffef4c 

}; 

Note: //Generate sine table from matlab// 

DSK6713_AIC23_Config config = { \ 

0x0017, /* DSK6713_AIC23_LEFTINVOL Left line input channel volume */ \ 
0x0017, /* 1 DSK6713AIC23RIGHTINVOL Right line input channel volume *A 
OxOOff, /* 2 DSK6713_AIC23_LEFTHPVOL Left channel headphone volume */ \ 
OxOOff, /* 3 DSK6713AIC23RIGHTHPVOL Right channel headphone volume 

*/\ 

0x00 10, /* 4 DSK67 1 3_AIC23_ANAPATH Analog audio path control */ \ 
0x0000, /* 5 DSK6713_AIC23_DIGPATH Digital audio path control */ \ 
0x0000, /* 6 DSK6713AIC23POWERDOWN Power down control */ \ 

0x0043 , /* 7 DSK67 1 3_AIC23_DIGIF Digital audio interface format */ \ 

0x008c, /* 8 DSK67 1 3AIC23S AMPLERATE Sample rate control */ 

\ 

0x0001 /*9DSK6713_AIC23_DIGACT Digital interface activation */ \ 

}; 
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DSK67 1 3_AIC23_CodecHandle hCodec; 

Int32InitWait=l; 

Int32 data; 

Int32 Logger[1024]; 

Int32 Loggerlndex =0; 

Int32 Avl =0; 

float in_buffer[ 100]; 



main() 

{ 



int i; 

LED=0x0; 

CSL_init(); 

// Filter Initialization 

for( i = ; i < 100; i++ ) in_buffer[i]=0.0; 

// Initialize codec 

hCodec = DSK6713_AIC23_openCodec(0, &config); 

IRQ_globalEnable(); 

IRQ_enable(IRQ_EVT_RINT 1 ); 

IRQ_enable(IRQ_EVT_XINT 1 ); 



void led(){ 

static int cc = 1 ; 

LED = cc; 

// To Shift Pattern 

if(cc== 0x03) cc = 0x05; 

else if (cc == 0x05) cc = 0x06; 

else if (cc == 0x06) cc = 0x03; 

else cc = 0x03; 
} 

setpll200M(){ 

} 

void read() { 

//int i =0; 
// float result; 

static int Index =0; 
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} 



data=MCBSP_read(DSK67 1 3_AIC23_DATAHANDLE); 
data += sinetable [Index]; 
Index+=2; 

if (Index >=48) Index %=48; 
Avl=l; 

Logger[LoggerIndex++] = data; 
if (Loggerlndex == 1024) 
Loggerlndex = 0; 



void write() { 

if(Avl=l){ 

MCBSP_write(DSK67 1 3AIC23DATAHANDLE, data); 
Avl =0; 

} 
else 

MCBSP_write(DSK67 1 3AIC23DATAHANDLE, data); 

} 



4. Enter the source code and save the file with main.c extension. 



'4 /Chipmax_6713/TMS320C6713_0 - C671x - Code Composer Studio - [main.c] 



^ File Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



» t£H| %d 



3 a a* «%« #*? 



tii ^i 



on|^%^%|^ 



|FIR6713M.pjt 



^||Debug 



~3\ * caaAl 4 *|^fc • £ 



SI oV 



ft 
{"'■> 

El 
n 

m 



a liiBBd 



J Files 
Ltl-Q GEL Piles 
Q -^ Projects 

!- 1 ^ FIR6713M.pjt (Deb 

C^| Dependent Project 
Pi Documents 

+ Cj dsp/bios config 

+ d Generated Files 
+ C3 Include 

CJ Libraries 
+ CfJ Source 

[^| boot.cmd 



-^r 



l#include <stdio.h> 

l#include "c6713dsk.h" 

l#include "master. h" 

l#include "aic23cfg.h" 

l#include "dsk(5713_aic23 .h" 

l#include <std.h> 

l#include <swi .h> 

l#include <log.h> 

l#include <c6x.h> 

l#include <csl .h> 

l#include <csl_mcbsp .h> 

L''* Length of sine wave table *s 
l#define SINE_TABLE_SIZE 48 

\ Delta 

Y''*floa t fil ter_ Coefff] ={8. OB ', 8. 80 ', O. 88, 8. 00, 8, 00, 8. 00, 8. 88, 8. 88, 

\8. 88, 8. 88, 8. 88, 8. 88, 8. 88, 8. 88, 8. 88, 8, 88, 8. 88, 8. 88, 8. 88, 8. 88, 8. 88, 8. 88, 8. 8i 

\8. 88, 8. 88, 8, 88, 8. 88, 8. 88, 8. 88, 1, 88, 8. 88}; *S 



Y's' Low Pass Filter 
V'* float filter_Coeff[] * 
\-l, 941786448798678e-882, 
\-6. 859448788578791e-882, 
1-6, 859448788578791e-882, 
\-l, 941786448798678&-882, 



■{2. 7152976341 71 146e-883, 2, 31581 9882942924e-884, -. 
-1, 72989821 8843226e- 882, 1, 773888738568675e-882, , 
8. 97831 '3256 '4 '4 826 Be - 882, 3. 81 4572949374625e - 881, 4. 
-6. 788468549771 738e- 882, 2. 11 3436751 136944e-883, -. 
2. 734 134 IB 7 e - 883, 1 . 28734 171 61 54 354e - 882, 7. 
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5. Right click on source, Select add files to project and Choose main.c file Saved 
before. 



_6713/TMS32QC6713_0 - C671x Code Co 



fy File Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 






D g£ B | & % B | ■ 



tj| £ ft "U ™ % u m *? tSi -5= | «!! &i | o n | a % % % | < 



|FIR6713M.pjt 



^|| Debug 



-T]| # m m *| ® «|^» • £ 



] ov aisiasd 



f Files 
+ CJ GEL files 
- £3 Projects 

B ^ FIR6713M.pjt (Debi 

I Dependent Project 

CD Documents 

+ C3 DSP/BIOS Config 

,+ C3 Generated Files 

+ CD Include 
Q Librarie s 

+ £2 — 

HbO( 




l#include <stdio.h> 

l#include "c&7 13dsk . h " 

#include "master. h" 

l#include "aic23cf cj . h " 

Itfinclude "dsk6713_aic23 .h" 

Itfinclude <std.h> 

l#include <swi .h> 

l#include <locj.h> 

l#include <c6x.h> 

l#include <csl .h> 

|#include <csl_mcbsp .h> 



f Lengtn of sine wave table *>'' 
lefine SIWE_TABLE_SIZE 48 






Beit a 
floa t. fii t.er_ Coefffj ={Ll. 88 ^ 8. Q8 r 8. 80 ', 8, Q8 r 8, QQ r 8. 80 ', tf„ £W„ £7. Q8„ 
8, 88^ 8. 80 ', tf. Z.7/.7, tf. tftf, tf . tftf, tf. Z.7/.7, tf, GG r 8. QQ r 8, QG r 8. QQ r 8. OQ r 8. 80 ', tf . £W, tf . ffi„ 8, GO, 8, GO, 
8, 80 '„ tf. £W, tf . tftf, tf. tftf, tf . tftf, tf. Z.7/.7, 7, tfff, 8. 88} ■ *.-** 

>' ZiW / 7 c?5 v 5 v Filter- 

'#fioa t fii t.ejr_ Coefffj ={2. 7152976341711 46e - 883 ', i\ 31 581 9882 942 924e - 684 r - 1 . 244493. 
1,941781 82^ -1.72981 ?26e-882^ 1 . 773881 ?5e-882 r 4.891495. 

6. 859441 91 e - 882 ', «?. 97831 3756 4 48266e - 882 r 3. 81 457294 9374625e - 881 ', *. 8198894. 

6. 859441 91 e - 882 ', - tf. 788468549771 738e - 882 ', i 7 , 11 34367511 36 944e - 885 ', #. AW #.95". 

1.941781 82 ', 2.734134 32e-883 r 1 . 287341 7161 54354e-882 r 7.2443649. 



7? £T?L?7^-Z 



_ V7V7 ■ J ^ . j= * 



^ 



1 



nect to each target . 
nect to each target . 



Sequence ID: 
Sequence ID: 



Error Code : 
Error Code : 



-1070 Error Class: 0x80000200 
-1070 Error Class: 0x80000200 



~~ 1 | | f\ Messiiyes /~~ 



ILd 



6. Add the other supporting .c files which configure the audio codec. 





Add Files to Project 


mm 




Lo 




3- m- 




ok in: |£3FIR6713M _J 'Sf 




(El 


Debug |Tl Q 
FIR6713M.CS_ IjO main.c 






IS 

e 

IS 

IS 


fcbap'" 








File 
File 








name: |"pll. c" "aic23.c" "aic23cfgLc.c" "dsk671 3_aic2 


Open | 




Cancel 




s of type: | c Source Files ( K .c; K .ccc) _^J 






Help 








d 
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7 Right click on source, Select add files to project and Choose boot.asm and save 
it. 



' /C671X AH_XDS510 Emulator/TMS320C671x_0 - C671x - Code Composer Studio - Not Connected - [main.c] j- ||n||Xl 



4p File Edit View Project Debug GEL Option Profile Ti 



& 



%a| 



~3h t*%^l#tf 



_^]| Debug 



Q i^L10® d 



+ _| GEL files 


B -^ Projects 


B ^ Test-enpt5.pit ( 




• Qj. Dependent Pro 




■■ CD Documents 




■■ Q D5P/BI05 Conf 




■■ Q] Generated File: 




■ [~~J Include 




C] Libraries 


B-^ Source 


! 


: ] aic23cfg_c 


! 


"| 3ic23cfg_c 


; 


S] boot.asm 


I 


"J dsi 6 i 3_a 


i 


,:\ dsk67i 3_a 


: 


] dsic7i.: a 


i 


S) main.c 


i 


1 pll.c 


v 


< , UN | > 





QM 



DISCONNECTED (UNKNOWN) 




dJ 



Look in: j £D Test-expt5 



S fji' 



!S) Test- expt5,C5_ 

B IW1M 



File name: 
Files of type: 



| boot 



Open 



Asm Source Files (".a") 



"3 



17 
18 
19 
20 
21: 



|2. 7 15297634 17 1146e-003}; 



-003, 
26e-0 
6e-00 
30e-0 
2e-00 



Iss High Pass Filter 
S*floa t til ter_ Coefff I = {3. 294316420702696s - 00 
mi. 3235* .■---. \ '.: : ■■. '.■■; ; ■ ie-004^ „ 



For Helpj press Fl 



JLn 1, Col ^ 



8 Right click on source, select add files to project and choose boot.cmd and save it. 



or/TMS320C671x_0 - C671x - Code Composer Studio - Not Connected - [main.c] |- 



4p File Edit View Project Debug 



i\ & y I * ii m h 



~T]| fa?ftt5 §» 



JTest-eHptS.p 



^ || Debug 



6V Q 11^0g£i 



+ _| GEL files 
B-^3 Projects 

B-^ Test-eKptS.pjt ( I 

: •• C7J Dependent Pro 
C~J Documents 
_| DSP/BI05 Conf 
• C~j Generated File:: 
I • C~j Include 
I -C~} Libraries 
B £3 Source 

~E\ aic23cfg_c 
^] aic23cfg_c: 
£] boot.asm 
\£\ dsk6713_a 
[1| dsk6713_a 
!■■ \£\ dsk6713_a 
£] main.c 
£] pll.c 
L?l boot.cmd ^ 

.<A mi | > 



Add Files to Project 



'JJ 



Look in: | tD Test-exptS 



Hf Stfi' 



||j^)Tesr_-expr_5,C5_ 



File name: boot 

Files of type: | Linker Command File p\cmd;Mcf] 



Open 



"31 



Help 



-003, 
26e-0 
6e-00 
30e-0 
2e-00 



■ 2.7 15297634 17 1146e-003}; 

l-'' ; ''' Hi$& Pass Fil ter 

\"* float filtoT_Cooff[j ={3. 294316423702696e-0l> 

i 1 . 32351 70 1 • 7$ / 4 go x 8e - 002 r 2, 635 9e. - 004, 



% »_■ DISCONNECTED (UNKNOWN) 



T 



For Help, press Fl 



|Ln lj Col J 
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9 Right click on DSP/BIOS Config and select add files project and add aic23.tcf file and 
save it. 



* /C671X AH_XDS51G EmuLatorJTMS320C671x_G C671x Code Composer Studio Not Connected [main.c] | - 



^ File Edit View Project Debug GEL Option Profile Tools D5P/BI05 Window Help 



J] & 



~T]|, 'ftVSfrS&lltf 



tj= +)= | H>= H>= | 



| Test-expt5.pjt 



^ || Debug 



£3 6V Q iB3[^E@d 



+ 


_| GEL Piles 




- 


£3 Projects 
- ^ Test-eHpt5.pjt ( 

l CD Dependent Pro 








Documents 




DSP/BIOS Conf 


1 aic23.tcf 




[+' EH Generated File 






I Include 






I Libraries 






- _J Source 






! | 


j*j| aic23cfg_c 






! ! 


jj| boot. asm 






| j 


1) dsk6713_a 






! ! 


J] dsk6713_a 






i i 


1) dsk6713 a 






! i 


lij main.c 






I ! 


1 pll.c 




i H boot.cmd 


< 


mi I > 



D/* 



< 



71 I f& fffi *±± j£ I Jl| Jig! I 



ft, 



lL^)Test-expt5.C5, 




] 



File name: 

Files of type: | Configuration File ( K .tcf) 



12. 7 15297634 17 1146e-003}; 

Y'-'y' High Pass Fil t&T 

\''*floa t fil t&r_ Co&fffJ ={3. 294316420702698& - 00 

\l,323l >-002 r 2„6356 s-Q84'J\ 



DISCONNECTED (UNKNOWN) 



~r 



For Help, press Fl 
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8 Right click on Libraries and select add files to project add csl6713.1ib from 
C:CCStudio v3.3/C6000/csl/lib/csl6713. 



* /C671X AHJCDS51Q Emulator7TMS320C671x_G - C671x - Code Composer Studio - Not Connected - [main. 



4^ Fil 



a) i-b 



Tifi:a'h^ | ft?i|i»iff 



(3 .+(= 

M:i= His 



jTest-enptS.pjt 



"3 [Debug" 



w 

?> 



{+} 



El 



ffl at o IS10I d 






Ei- ■■□ GEL files 
B- S Projects 

B- jS Test-eKpt5.pjt ( 

' -£~} Dependent Pro 

L"J Documents 
B- £3 DSP/BIOS Conf 

■ aic23.tcf 
[+• \^~} Generated File: 

: Q Include 

B ^ Libraries 
B-Q Source 

| [3 aic23cfg_c: 

=£] boot. asm 
£j dsk6713_a 
J] dsk6713_a 
£j dsk6713_a 
;■■ [J] main.c 

2 pii.c 

fil boot.cmd 

> 



< 



-Dj>r 



■ DISCONNECTED (UNKNOWN) 



LJJ 



"1 1 ft* iiii ±±± 



Add Files to Project 



A * I 



@S 



Look in: Q lib 



"3 * S8fr 



|^csl6413 

Hcsl6413e 

[H)csl6414 

[ffijcsl6414e 

Hcsl6415 

Hcsl6415e 

m 



j(csl6416 

jjcsl6416e 

jjcsl641S 

j|csl6413e 

jjcsl6701 

f)csl6701e 



l^csl6711 

Hcsl6711C 

Hcsl6711Ce 

[EJcsl6711e 

Hcsl6712 

Hcsl6712C 




File name: Icsl6713 

Files of type: | Object and Library Files [".o^T] 



Open 



"3 



Cancel 



Help 



003, 

f26e-0 

6e-00 

30e-0 

2e-00 



17 
18 
19 
20 
21 



|2. 7 15297634 17 1146e-003}; 

I,'-;-'' High Pass Fil t.&jr 

>s'* -float. filt.&r_Co&ff[j ={3,294316428702696&-0Q 
MI, 323547807544983&-002^2, 63589698604891 9&-804, 



For Help, press Fl 



|Ln 1, Col ^ 
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9 Go to Project- -> Build option -> Preprocessor^ set CHIP6713 as shown below. 



*■ /C671X AH_XDS510 EmulatorMMS320C671x_0 - C671x - Code Composer Studio - Not Connected - [csl_chiphal.h] 



y File Edit Vev Project Debug GEL Option Profile Tools DSP/BIOS Window Help 



I g£ B I * iH e | ' 



| Test-Code-Del. pit _▼][ Debug 



-gi^niii^l^^ +%> 



to 

m 

m 
n 
m 
m 



a isiBSd 



$ Files 

+ _| GEL files 

- £3 Projects 

- ^ Test-Code-Del.p 

!■■ CH Dependent Pro 
i--Pl Documents 
El ^ DSP/BIOS Conf 

i--EJl aic23.tcF 
El £3 Generated File: 
[J] aic23cfg.se, 
[1] aic23cfg_c 
LJ Include 
:-: ^ Libraries 

■■\£\ csl6713.lib 
;-: £3 Source 

^ aic23.c 
£] aic23cfg_c: 
:tn boot, asm 
I) dsk6713_a 
1 dsk6713_a 
1 dsk6713_a 
1] main.c 
Id pll.c 
■ U^\ boot.cmd 



T]| ft ft 'fa w % *4 a *? t!a +5i I «£a tf | O H | /i r * & % | < 



251 
25 
251 
25' 
261 
26 
26 
26 
26 
26 
261 
26 
261 
26 
271 
27 
27 
27 
27 
27! 
271 
27 
271 
27' 
231 



! General Compiler | Linker j DspBiosBuilder | Link Order | 



■g -pd • '22F. -fr"$(F :-] 3n!'.Dt-t ig d"_L t 3UG ' nv iV'.'n 



Category: 


Preprocessor 


Basic 

Advanced 

Feedback 

Files 

Assembly 

Parser 

Diagnostics 


Include Search Path (-i): | 


Pre-Define Symbol (-d): |CHIP_6713J 


Undefine Symbol (-u): \ 


Preprocessing: JNone _^J 
|~~ Continue ■■ it' C mpilation (-ppa) 
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sre XXXX is chip number, i.e. 6; 



cS,c9,clO,cll,cl2,cl3,cl4,cl5,: 



.- 



1 fatal error detected in the compilation of "main.c 
Compilation terminated. 

>> Compilation failure 



| |\ Messages^ Build f~ 



10. a) Go to Project to Compile. 

b) Go to Project to Build. 

c) Go to Project to Rebuild All. 



6713/TMS320C6713 0-C671x-C 



uJiL.ujjjmummi.il 



► File Edit View 



I Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



% & y | J 



|FIRG713M.pjt 



^ 6V Q I 



Open... 

Use External Makefile... 

Export to Makefile... 

Add Files to Project... 

Save 

Close 




Configurations... 
Build Options,., 

File Specific Options... 



Project Dependencies. . . 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 



File View I /4 Bookmarks 



J | tf ft "fo <M L 


s* tk | m *? tsi -si 


m m s\ * m\^% 


• £ 


■ #B|Q 



include <stdio.h> 
include "c6713dsk.h" 
#include "master. h" 
#include "aic23cf g . h " 
#include "dsk6713_aic23 .h" 
^include <std.h> 
#include <swi .h> 
#include <log.h> 
#include <c6x.h> 
#include <csl .h> 
#include <csl_mcbsp .h> 

'* Length of sine wave table *,--' 
tfdefine SIWE_TABLE_SIZE 46 

' Delta 

■■floa t fil ter_ Coefffj ={8. 88, 8. 80 ', 8. t 
8. 88, 0. 00 r 0. 88, 0. 00 r 0. 00 ', 0. 88, 0. 00 ', 0. t 
0, 00 ', 8. 00, 8, 88, 8. 80, 8, 88, 8, 88, 1, 88, 8. , 

''.••' Low Pass Filter 

-'* 'float filter_Coeff[j ={2.7152976341, 
-1. 941786- le-882, -1. 72989821 8, 

-6. 859448, te-882, 8. 9783132564- 

-6. 859448, e-882, -€. 788468549, 

-1, 941786-, ?e-882, 2. 7341341665, 
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d 



< 



D File View /4 Bookmarks 



\s* float fllt&r_Coeff[] ={2, 715297634171 146o-883^ 2. 31581980294292 



78&-882 t 
I - 6. 859448788578791& - 882 r 8, 9783132. 
\-6. 8594 48788578 79 1&- 882 * 
\-l. 941 786 4 487986 78& - 882 ', , 



'374625 
7113694 
154354 



[Linking. . . ] "C:\CCStudio_v3. l\C6000\cgtools\bin\cl6x^ 
<Linking> 

Build Complete, 

Errors, Warnings, Remarks. 



*;/////£% It will shows 
warnings and 
errors if any 



| | |\ Messages )\ Build f~ 
%g ■a|HALTED 



1 1 . Put a break point at Loggerlndex = 0; 

12. Go to file and load program and load ".out" file into the board. 



ffiW.ULI J J..UUfJ.ULl t BJ.lll«J.IJJJ.IJJ.UJlllll.— M 



Edit View Project Debug fciEL Option Profile Tools PB>:_ DbP/BIOb Window Help 



Close 
Save All 





J | #4 Ml *%* M 


% t.f B | m *? 


a ^ 


44 = 


*1 


|00|^3 


^l^n 


tt « * 1^ % 


• *| 










I £fr lSj 


Hia 





Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 



Load GEL... 

Workspace 
File I/O... 



Difference betwee 

Merge Fibs.. . 



Print.. 



Recent Source Files 
Recent Workspaces 

Recent Program Files 
Recent Symbols 
Recent GEL Files 



Launch betup 



#include <stdio.h> 

l#include "c6713dsk.h" 

l#include "master. h" 

l#include "aic:23c:f g . h " 

l#include "dsk67 13_aic23 . h " 

l#include <std .h> 

#include <swi .h> 

;lude <log.h> 

#include <c6x.h> 

|#include <csl .h> 

|#include <csl_mcbsp.h> 

'''* L&nfft-A oF sine wave taitl& *=-•'' 
|#define SIHE_TABLE_SIZE 4S 

:-'' Belts 

'*Flcra t. Fll t.&jr_ CosFFfJ ={G. GO, G. GO, G. GO, O, GO, G. GO, G, G6 
iO. GO, O, GO, G. GO, G. GO, O. GO, G. GO, G. GO, G. GO, O. GO, G. GO, G, G5 
\o. OG, O. GO, G. GO, G, GO, O. GO, G. GO, 1. GO, G. GO J,- *,--' 

.'•','•' Law Pass Fj'lt&JT 

,--'*Flaat. F2lt&jr_Ca&FffJ ={2.71. ?-OG3, 2. 



-6. 93944 7 91e-992 

-6. 9394 ■ 791&-OG2 

-1.94171 ?&-OG2 



315 
-1.729G& 92,1.775 

8, 979315 'Be - OG2, 3. Gl 45 

-6.7884b 92,2.115 

2. 734134 3, 1. 2G75 



Load Program 



Look in: | ED Debug 



~z\ <= s & nnn- 



File name: Iaic23.< 

Files of type: I K .out 



Open 



"3 



Cancel 



Help 



T7^7lJ 
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13. Go to Debug and click on run to ran the program. 



1 GEL Option Profile lools PBC D5P/BIOS Wi 




UUUUMUH UUUUUb^ 

nnnnvAnn n 1 DntnoA 



14. To see the Graph, go to View and select time/frequency in the Graph, and give the 
correct Start address provided in the program, Display data can be taken as per user. 



E Graph Property Dialog 




*] 


Display Type 


| Single Time 


Ml^ 


Graph Title 


Graphical Display 




Start Address 


Logger 




Acquisition Buffer Size 


1024 




Index Increment 


1 




Display Data Size 


1024 




DSP Data Type 


32-bit signed integer 




Q -value 







Sampling Rate (Hz) 


16000 




Plot Data From 


Left to Right 




Left-shifted Data Display 


Yes 




Autoscale 


On 




DCValue 





— ' 


Axes Display 


On 




Time Display Unit 


s 




Status Bar Display 


On 


T | 








Help | 


| OK Cancel | 




115 



15. To view the Spectrum, go to View and select time/frequency in the Graph, and give 
the correct Start address provided in the program, Display data can be taken as per user. 







x] 


OHHS 


FFT Magnitude 


d* 


Graph Title 


Graphical Display 


Signal Type 


Real 


Start Address 


Logger 


Acquisition Buffer Size 


1024 


Index Increment 


1 


FFT Framesize 


1024 


FFT Order 


10 


FFT Windowing Function 


Rectangle 


Display Peak and Hold 


Off 


DSP Data Type 


32-bit signed integer 


— i 


Q-value 





Sampling Rate [Hz] 


16000 


Plot Data From 


Left to Right 


Left-shifted Data Display 


Yes 


Autoscale 


On 


d 








| OK | Cancel | 


Help | 





3.3e+7 
3.1e+7- 
2.8e+7- 
2.5e+7- 
2.2e+7- 
1 .9e+7- 
1 .7e+7- 
1 .4e+7- 



5.6e+6- 
2.8e+6- 



vW 



wj 



Vy^^^y^ y^,^r^4r/^ ,v-; '' V v ^^?^^^^ 



667 1333 2000 2667 3333 4000 4667 5333 6000 



(7984,38,660702) 



|FFT Mag 



Lin Auto 5cale Rectangular 
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Experiment 6: Noise Removal 



Procedure to create new Project: 

1. To create project, Go to Project and Select New. 



3 /ChipMax_671 37TMS320C671 3_0 - C671x - Code Composer 



File Edit View 



£ g£ 



r 



£3 &r 



{? 



jf Files 
B-C] Gl 



Debug GEL Option Profile Tools PBC DSP/ 



Open.,. 

Use External Makefile. . . 

Export to Makefile... 

Add Files to Project... 

Save 



Sttti 



Source Control 



Compile File 
Build 

Rebuild All 
Stop Build 
Build Clean 



Ctrl+F7 
F7 




2. Give project name and click on finish. 



^^^^^^ 








Project Name: | 








Location: | Q \CCS tudio_v3. 1 \MyPrcjects\ 


_l 






Pi ujeut Type. | E xecutable (. out) 


jd 






Target |tMS320C67>C< 


zl 






< Back Finsh 


Cancel 


Help 






^ 



(Notej^ Location must be c:\CCStudio_v3.3\MyProjects). 

3. Click on File ' ^ New ' ^ Source File, to write the Source Code. 



Edit View Project Debug GEL Option Profile Tools PBC DSP/BIO^ 



Open,.. 


Ctrl+O 








Ctrl+S 


5ave As,.. 




5ave All 




Load Program.,. 


Ctrl+L 


Reload Program 


Ctrl+5hift+L 


Load Symbols 


► 


Reload Symbols 


► 


Unload Symbols 


► 
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AIM: Noise removal in a given mixed signal. 

Noise removal: In the real time systems every signal will get corrupted by the noise. To 
analyze the noise we need to add the noise, a noise will be generated and added to the 
signal. A noise of 3KHz is added to a tone of 400KHz then the noise is removed by using 
an high pass filter. 

Program: 

#include <stdio.h> 
#include "c6713dsk.h" 
#include "master.h" 
#include "aic23cfg.h" 
#include "dsk6713_aic23.h" 

#define FiltLen 47 

#define SINE_TABLE_SIZE 48 

float PrNoise 

[600]={0.3509,0.1 118,0.3631,0.2673,0.2826,0.1909,0.2901,0.1378,0.1270,0.1798, 

0.2925,0.3730,0.1307,0.1380,0.2859,0.1738,0.3278,0.0355,0.3177,0.2905, 

0.1473,0.2380,0.2112,0.2662,0.3228,0.0264,0.2991,0.3027,0.1220,0.1676, 

0.2360,0.2942,0.3058,0.3013,0.1423,0.1478,0.3479,0.2017,0.3278,0.1164, 

0.2747,0.2752,0.4239,0.0666,0.2130,0.2124,0.0124,0.3306,0.3134,0.1420, 

0.2837,0.0446,0.2438,0.1069,0.1902,0.3497,0.0962,0.1404,0.2046,0.1419, 

0.3231,0.3962,0.0549,0.0875,0.2742,0.1418,0.3797,0.3756,0.2441,0.3271, 

0. 1456,0.2332,0.3789,0.2250,0.3798,0.3 1 53,0.2295,0.2950,0.3924,0.0982, 

0. 1493,0.3649,0. 1 1 59,0.2095,0.3088,0.3656,0.2539,0.0606,0.0548,0.2847, 

0.2610,0.4005,0.2985,0.2459,0.1789,0.0824,0.3774,0.2256,0.0309,0.2949, 

0.2133,0.0136,0.3288,0.0923,0.0491,0.0178,0.1205,0.1107,0.3042,0.0761, 

0.1512,0.1302,0.3434,0.1384,0.2454,0.1873,0.2804,0.0407,0.4164,0.3858, 

0.0678,0.2275,0. 1430,0.2304,0.0500,0. 1 793,0. 1 887,0.2067,0. 1 776,0.0374, 

0.2449,0.2093,0.0772,0.1982,0.1787,0.2021,0.2916,0.2557,0.3304,0.0184, 

0.0892,0.1823,0.3762,0.2984,0.1576,0.1120,0.0088,0.3026,0.2274,0.1223, 

0.2151,0.1131,0.3248,0.1441,0.3475,-0.0471,0.2749,0.1925,0.1708,0.2431, 

0.2975,0.0634,0.3891,0.0372,0.1486,0.0943,0.3835,0.0355,0.2320,0.2466, 

0.2428,0.3175,-0.0105,0.2891,0.1764,0.2621,0.0814,0.2239,0.3074,0.4196, 

0.2065,0.0613,0.0321,0.2495,0.3044,0.2513,0.1193,0.3635,0.2031,0.3801, 

0.0303,0.0733,0.3001,0.3507,0.2985,0.1186,0.0656,0.3612,0.3397,0.1294, 

0.1102,0.1194,0.3025,0.2892,0.1845,0.1956,0.1073,0.2981,0.3682,0.2311, 

0.2244,0.2099,0.0990,0.2220,0.2894,0.2813,0.2316,0.0653,0.2872,0.1048, 

0.1335,0.1639,0.1335,0.2752,0.0262,0.2958,0.2226,0.2916,0.1142,0.2017, 

0.3252,0.2093,0.1280,0.4335,0.0627,0.1850,0.4388,0.1821,0.1451,0.0544, 

0.3462,0.0154,0.0822,0.3031,0.1478,0.2130,0.2230,0.2897,0.0397,0.2436, 

0.0428,0.3551,0.1458,0.3382,0.0957,0.2303,0.3804,0.0262,0.0356,0.0130, 

0.2607,0.2102,-0.0020,0.2581,0.1933,-0.0047,0.0244,0.0922,0.3805,0.1092, 

0.1670,0.1447,0.0477,0.3077,0.2124,0.2124,0.1879,0.1623,0.1219,0.2904, 
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0.0953,0.1132,0.0502,0.4002,0.3765,0.0096,0.1440,0.2291,0.1027,0.3114, 

0.2580,0.2089,0. 1434,0.3 168,0.3503,0.255 1,0. 1064,0.3 127,0.0588,0. 1926, 

0.3867,0.1418,0.1538,0.4241,0.0659,0.1438,0.2286,0.0508,0.2072,0.2740, 

0.1688,0.2416,0.0827,0.2314,0.3496,0.1934,0.3886,0.0894,0.0704,0.1493, 

0.1843,0.1859,0.1252,0.2559,0.0035,0.1217,0.0717,0.0912,0.0251,0.2405, 

0.1436,0.1074,0.0493,0.1552,0.2456,0.3565,-0.0167,0.3025,-0.0187,0.2772, 

0.2816,0.3582,0.0750,0.2422,0.2169,0.1210,0.2634,0.2424,0.0600,0.1631, 

0.3293,0.1383,0.2371,0.1551,0.0338,0.1593,0.3720,0.1812,0.0607,0.1999, 

0.1206,0.2411,0.1635,0.2727,0.2958,0.0758,0.0701,0.3565,0.2880,0.0363, 

0.3726,0.2002,0.2003,0.2127,0.2768,0.3146,0.1400,0.0291,0.2936,0.1341, 

0.3375,0.1544,0.3321,0.0716,0.0532,0.3473,0.0176,0.2671,0.2772,0.1617, 

0. 1264,0.3688,0. 1475,0. 1768,0.0764,0. 1556,0. 1539,0.3684,0.0979,0. 1012, 

0.1261,0.2401,0.0153,0.3234,0.0373,0.2052,0.0465,0.3405,0.3056,0.0901, 

0.2585,0.1746,0.2905,0.1795,0.3366,0.1744,0.1618,0.2372,0.1421,0.3937, 

0. 1 927,0.0760,0. 1 248,0.2367,0. 1 1 59,0.043 1 ,0.3350,0.2452,0. 1 800,0. 1 234, 

0.1133,0.2164,0.1742,0.2783,0.0053,0.3180,0.2518,0.0386,0.3270,0.0609, 

0.2273,0.2060,0.0477,0.0029,0.3 1 82,0.32 1 8,0. 1 980,0.0483,0.2532,0.0704, 

0.2688,0.1805,0.0634,0.3304,0.2816,0.3470,0.0396,0.0822,0.3077,0.2812, 

0.2906,0.1659,0.1466,0.2278,0.3560,0.3095,0.2671,0.1798,0.3339,0.1612, 

0.1967,0.2755,0.2225,0.2483,0.3013,0.2941,0.0201,0.0807,0.1395,0.2114, 

0.0911,0.2137,0.1130,0.1435,0.2817,0.0310,0.0678,0.3107,0.1429,0.0814, 

0.2429,0.3712,0.1587,0.2622,0.2169,0.1503,0.1834,0.3897,0.0846,0.1703, 

0.2341,0.2631,0.3412,0.1344,0.0903,0.1993,0.0939,0.2847,0.1697,0.2312, 

0.2761,0.1753,0.1190,0.0868,0.3864,0.0813,0.3023,0.2145,0.2764,0.3148, 

0.0086,0. 1 329,0.3952,0.2684,0.3843,0.2 1 8 1 ,0.0226,0.268 1 ,0. 1 080,0.2249, 

0.2832,0.2395,0.2193,0.1590,0.0663,0.1474,-0.0163,0.2260,0.3727,0.0303, 

0.0235,0.3627,0.1062,0.0084,0.3246,0.0568,0.3485,0.0407,0.1353,0.2338, 

0.2315,0.2545,0.2482,0.1510,0.0363,0.2882,0.0925,0.2603,0.3444,0.0445, 

0.4034,0.0162,0.3319,0.3212,0.0066,0.2010,0.3604,0.1120,0.3155,0.2390, 

0.3231,0.3301,0.0269,0.1138,0.2520,0.1875,0.3778,0.2939,0.2133,0.1170, 

0.0977,0.3524,0. 1 664,0.32 1 1 ,0.0252,0.336 1 ,0.0023,0.35 1 3,0.3686,0.2 151, 

0.0581,0.0777,0.1664,0.3211,0.0252,0.3361,0.0023,0.3513,0.3686,0.2151}; 

Note: // Generate the Noise from the Noise table from the matlab using rand() 

function// 

float RemNoise[47]={ -2.181030629062908e-002,8.917428211564256e- 

003,1. 182863 130963947e-002, 1.428453375748106e-002,1.403350814552515 e , 

002,1.039856067649067e-002,5.093274200275827e-003, 1.421581262318843e- 

003,2.5 15430936577368e-003,9.3630532 131 1 1 1 26e-003, 1.9491 52 1 12446623e-002, 

2.722897698684972e-002,2.5799432377005 17e-002, 1 .050645 1 9 1 333675e-002, 

1.8321 63998604397e-002, -5.406 1 941 8 1 033640e-002,-8.5 148697235846 1 6e-002, 

9.887618777194764e-002,-8.633153232820758e-002,-4.651515024517261e 

002, 1 .2 1 70 1 1 6 1 0629542e-002,7.394838432088444e-002, 1 .2068026 1 6409422e-00 1 , 

1 .380624377729094e-00 1 , 1 .2068026 1 6409422e-00 1 ,7.394838432088444e 

002, 1 .2 1 70 1 1 6 1 0629542e-002, -4.65 1 5 1 50245 1 726 1 e-002,-8.633 1 53232820758e-002, 

9.887618777194764e-002,-8.514869723584616e-002,-5.406194181033640e-002,- 
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1.832163998604397e-002,1.050645191333675e-002,2.579943237700517e-002, 
2.722897698684972e-002, 1 .949 1 52 1 12446623e-002,9.3630532 1 3 1 1 1 126 e 
003,2.5 15430936577368e-003, 1.421581262318843e-003,5.093274200275827 e 
003,1.039856067649067e-002,1.403350814552515e-002, 1.428453375748106e- 
002, 1 . 1 82863 1 30963947e-002,8.9 1 74282 1 1 564256e-003,-2. 1 8 1 030629062908e-002 } ; 

DSK6713_AIC23_Config config = { \ 

0x0017, /* DSK6713_AIC23_LEFTINVOL Left line input channel volume */ \ 
0x0017, /* 1 DSK6713_AIC23_RIGHTINVOL Right line input channel volume *A 
0x00d8, /* 2 DSK6713AIC23LEFTHPVOL Left channel headphone volume */ \ 
0x00d8, /* 3 DSK6713_AIC23_RIGHTHPVOL Right channel headphone volume */ \ 
0x00 10, /* 4 DSK67 1 3AIC23ANAPATH Analog audio path control */ \ 
0x0000, /* 5 DSK6713_AIC23_DIGPATH Digital audio path control */ \ 
0x0000, /* 6 DSK6713AIC23POWERDOWN Power down control */ \ 

0x0043 , /* 7 DSK67 1 3_AIC23_DIGIF Digital audio interface format */ \ 
0x008c, /* 8 DSK67 1 3AIC23S AMPLERATE Sample rate control */ \ 

0x0001 /*9DSK6713_AIC23_DIGACT Digital interface activation */ \ 

}; 

DSK67 1 3_AIC23_CodecHandle hCodec; 

Int32InitWait=l; 

Int32 data; 

float in_buffer[ 100]; 
float InSigBuffer[600]; 
float NoiseBuffer [600]; 
float CorrSigBuffer[600]; 
float RecovSigBuffer[600]; 
int Loglndex =0; 

main() { 

int i; 

LED=0x0; 

CSL_init(); 

for( i = ; i < 100; i++ ) in_buffer[i]=0.0; 

hCodec = DSK6713_AIC23_openCodec(0, &config); 

IRQ_globalEnable(); 

IRQ_enable(IRQ_EVT_RINT 1 ); 

IRQ_enable(IRQ_EVT_XINT 1 ); 
} 

void led() { 

static int cc = 1 ; 

LED = cc; 

if(cc== 0x03) cc = 0x05; 

else if (cc == 0x05) cc = 0x06; 
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else if (cc == 0x06) cc = 0x03; 
else cc = 0x03; 

} 
setpll200M(){ 

} 

void read() { 

int i = 0; 

float result; 

float InData; 

float InNoise; 

float InCorrup; 

static int NoiseIndex=0; 

data=MCBSP_read(DSK67 1 3_AIC23_DATAHANDLE); 

if(data>=32768) data= data|0xffffD000; 

InData = (Int32) (data)/ 16; 

InData=data; 

InNoise = PrNoise[NoiseIndex++]*4096 - 1024; 

if(NoiseIndex==600) 

Noiselndex = 0; 
InCorrup = InData + InNoise; 

for( i = 0; i <= (FiltLen-2); i++) in_buffer[i] = in_buffer[i+l]; 
in_buffer[FiltLen- 1 ]=InCorrup; 
result = 0; 

for( i = ; i <= (FiltLen-1); i++ ) result += RemNoise[i] * in_buffer[i]; 
data = (Int32)(result*512); 
InSigBuffer[LogIndex]= InData; 
NoiseBuffer[LogIndex]= InNoise ; 
CorrSigBuffer[LogIndex]= InCorrup; 
RecovSigBuffer[LogIndex] = (float)data; 
LogIndex++; 
if(LogIndex==600) 

Loglndex =0; 

} 

void write() { 

if(InitWait<1000){ 
InitWait++; 

MCBSP_write(DSK67 1 3AIC23DATAHANDLE, 0); 
MCBSP_write(DSK6713_AIC23_DATAHANDLE, 0); 



} 
else{ 



MCBSP_write(DSK67 1 3AIC23DATAHANDLE, data); 
MCBSP_write(DSK6713_AIC23_DATAHANDLE,data); 
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Output: 

Recovered signal can be seen as shown below. 




4. Enter the source code and save the file with main.c extension. 



'< /Chipmax_6713/TMS32QC6713_0 - C671x - Code Composer Studio - [main.c] 



^ File Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



i t£Q| %d 



3 a a.* ««%« #*? 



tii ^i 



on|^%%%|^ 



|FIR6713M.pjt 



^||Debug 



"31 1 * caaAl 4 *|^fc • £ 



SI oV 



ft 
{■> 

m 
n 



a liiggd 



f Files 
Ltl-Q GEL Piles 
B-£a Projects 

-^FIR6713M.pjt(Debi 

Q Dependent Project 

r~l Documents 
+ |~J DSP/BIOS Config 
+ d Generated Files 
E-Q Include 

CJ Libraries 
+ CJ Source 

[^?1 boot.crnd 



-oj>r 



l#include <stdio.h> 

l#include "c6713dsk.h" 

l#include "master. h" 

l#include "aic23cfg.h" 

l#include "dsk(5713_aic23 .h" 

l#include <std.h> 

l#include <swi .h> 

l#include <log.h> 

l#include <c6x.h> 

l#include <csl .h> 

#include <csl_mcbsp .h> 

L''* Length of sine wave table *v 
l#define SINE_TABLE_SIZE 48 

\ Delta 

\<''*floc? t fil ter_ Coefff] ={0. 00, O. 00, O. 00, O. 00, 0, 00, O. 00, O. 00, O. 00, 

\0. 00, O. 00, O. 00, O. 00, O. 00, O. 00, O. 00, O. 00, O. 00, O. 00, O. 00, O. 00, O. 00, O, 00, O. Oi 

\o. 00, O. 00, O, 00, O. 00, O. 00, O, 00, 1, 00, O. 00}; *S 



■J_r 



Y's' Low Pass filter 
V'* float filter_Coeff[] * 
\-l, 941796446799678e-002, 
1-6, 059440700570791e-002, 
1-6, 059440700570791e-002, 
\-l, 941 7064 40 7906 78e- 002 r 
7 it;? a -7/,- -i.4 ivii .*£„ _ nn ? ? 



■{2. 71 52976341 71 146e- 003, 2, 31581 9802942924e-004, -. 
-1, 72909821 38432260- 002, 1, 773008730568675e-002, , 
8. 97031 3256 4 48266e - 002, 3, 01 4572949374625e - 001, 4. 
-6. 788468549771 730e- 002, 2, 11 3436751 136944e-003, -. 
2. 73413416 7 e - 003, 1 . 20734 171 61 54 354e - 002, 7. 
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5. Right click on source, Select add files to project and Choose main.c file Saved before. 



6713/TMS320C6713 0-C671x 



^ File Edit View Project Debug GEL Option Profile Tools PBC D5P/BI05 Window Help 

1 & Q | 



^m 



|FIRG713M.pjt 



^|| Debug 



-T]l#aa^| ® &\^% +%> 



m at o isimsgii -%m 



f Files 

+ Cj GEL files 

- £3 Projects 

B -£§ FIR6713M.pjt (Debi 

I CJ Dependent Project 
| Documents 

+ Q DSP/BIOS Config 

Fy • O Generated Files 

+ CJ Include 
P~l Libraries 

+ CJ 

: dbo ( 



tj| iff ft •*» « % ?j a *? ts -si | H Si j@ | o n | > a % % I < 



l#include <stdio.h> 

l#include "c6713dsk.h" 

l#include "master. h" 

l#include "aic23cfcj . h " 

Itfinclude "dsk6713_aic23 .h" 

l#include <std.h> 

l#include <swi .h> 

l#include <log.h> 

Itfinclude <c(5x.h> 

l#include <csl .h> 

|#include <csl_mcbsp .h> 



i/ Allow Docking 
Hide 
Float In Main Window 



< 



Length of sine wave table *,-'' 
lefine SINE_TABLE_SIZE 48 



Delta 

bfloa t fil ter_ Coefff] ={8. 88, O. 88, 8, 88, 8, 88, 8. 88, 8. 88, 8. 88, 8. 88, 
\0. 88, 8, 88, 8. 88, 8. 88, 8. 88, 8. 88, O. 88, O. 88, O. 88, O. 88, 8. 88, 8. 88, 8. 88, 8. 88, 8, , 
\8, 88, 8. 88, 8. 88, 8. 88, 8. 88, 8. 88, 1, 88, 8. 88} ,■ */ 

■-''/' Low Pass Filter 

■<"*floa t fil ter_ Coefff] ={2. 7152976341711 46e - 883, 2. 315819882942924e - 884, 
- 1 , 94178 ~78e - 882, - 1 , 72989821 884 322 6e - 882, 1 . 7738887385686 75e - 882, 

- 6, 859448788578791& - 882, 8. 9783 J ' 32 56 4 4 826 6e - 882, 3, 81 4 57294 9374 625e - 881, 
- 6, 859448788578791e - 882, - 6. 788468549771 738e - 882, 2. 1134367511 36 944e - 883, 
- 1 , 94176 '78e - 882, 2. 7341 341 66585232e - 883, 1 , 287341 71 61 54354e - 882, 



-1. 244493. 

4. 891495. 
4, 8198894. 

4. 891495. 
7, 2443649. 



71 ^?a7£-3.4 



£s>-nni ? • 



7> 






nect to each target . 
nect to each target . 



Sequence ID: 
Sequence ID: 



Error Code : 
Error Code : 



-1070 Error Class: 0x80000200 
-1070 Error Class: 0x80000200 



TTT1V 



lid 



6. Add the other supporting .c /.asm/.cmd/.tcf files which configure the audio codec. 
Follow the same procedure which you have done in experiment 4 and 5. 

7. a) Go to Project to Compile. 

b) Go to Project to Build. 

c) Go to Project to Rebuild All. 




|FIRG713M.pjt 
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&\ o i 








g Files 








- Q Projec 


CJ* 


?} 


£ 


{? 


+ 


L 


a-c 


-r\ 


s-E 


rt 


_ 


^ 




2 






m 


Fl 






n 


1 


m 





Use External Makefile, 
Export to Makefile.,. 
Add Files to Project.,. 



Source Control ► 


-'"Compile Filers. 


Ctrl+F7 


Build ) 


F7 


-^J}ebuildAII_____-^ 




Build Clean 





Configurations.. 
Build Options.,. 
File Specific Options.., 

Project Dependencies. . . 
Show Project Dependencies 
Show File Dependencies 
Scan All File Dependencies 
Recent Project Files 



TT^ 



aa | * &\^% 



^include <stdio.h> 
#include "c6713dsk.h" 
#include "master. h" 
#include "aic23cfg .h" 
#include "dsk6713_aic23 .h" 
include <std.h> 
^include <swi .h> 
#include <log.h> 
#include <c6x.h> 
include <csl .h> 
#include <csl_mcbsp .h> 



:--" Pelt a 

'*flaa t fil ter_ Coefff] ={0, OO, O, OO. 
O. OO, O, OO, O. OO, O, OO, O. OO, O. OO, O. OO. 
O. OO, O, OO, O. OO, O. OO, O. OO, O. OO, 1. OO. 

' low Pass Filter 
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d 



< 



D File View /4 Bookmarks 



\s* float fllt&r_Coeff[] ={2, 715297634171 146&-883, 2. 31581980294292 
- 1 . 941 786 4 487986 78& - 882, - 1 . 72989821 8843226e - 882 ', 1 . 7738887385686/ 



91&-882, 3,9783132. 
\-6. 8594 48788578 79 1&- 882, 
\-l. 941 786 4 487986 78& - 882 ', , 



[Linking. . . ] "C:\CCStudio_v3. l\C6000\cgtools\bin\cl6x^ 
<Linking> 

Build Complete, 

Errors, Warnings, Remarks. 



? >S9 It will shows 
warnings and 
errors if any 



"374625 
7113694 
154354 



| | |\ Messages y\ Build y^ 
%g ■a|HALTED 



8. Put a break point at Loglndex = in the source code. 

9. Go to file and load program and load ".out" file into the board. 



6713ATMS320C6713 0-C671* 



Edit View Project Debug GEL Option Profile Tools PBC DSP/BIOS Window Help 



New 
Open,,. 
Close 
Save 
Save As,, 



Ctrl+O 
Ctrl+S 



^^^^^^ilQ 





IBEfflJJ.I.IJ.1..— 


Reload Program 
Load Symbols 
Reload Symbols 
Unload Symbols 


Ctrl+5hift+L 

► 

► 


Load GEL,,, 
Data 

Workspace 
File I/O,,, 


► 


Difference between files, . . 
Merge Files,.. 


Print,,, 


Ctrl+P 



Recent Source Files 
Recent Workspaces 
Recent Program Files 
Recent Symbols 
Recent GEL Files 



Launch Setup 



T* 



TJ| £ ft fl k w % U m iff <& +% | «5i &l I El H I A 



"311^ Bit* I €5 #|>^ •£ 



l#include <stdio.h> 

l#include "c6713dsk.h" 

l#include "master. h" 

l#include "aic23cfcf .h" 

Itfinclude "dsk6713_aic23 .h" 

l#include <std.h> 

l#include <swi .h> 

l#include <log.h> 

l#include <c6x.h> 

l#include <csl .h> 

l#include <csl_mcbsp .h> 

|--"* Length of sine wave taMe #S 
Itfdefine SINE_TABLE_SIZE 48 

I ; Delta 

s*floa t £11 ter_ CoefffJ ={8, 80 ', O. 00, 0, 00 ', tf , tftf, tf. 00, O, 06 
\0. 88, Q. 88, Q. 00 ', A tftf, A tftf, tf, 88, 8. 88, 8, 80 ', tf , tfZ7, tf. ffi7„ tf, tf£ 
L?- tftf, tf. 88, 8, 88, 8. 88, 8. 88, 8, 88, 1. 88, 8, 88/; *,■" 

\' r '--'' Low Pass Filter- 

\-'* float fllter_Coeff[J ={2, 71 52976341 711 46e- 883, 2,315 

\-l, 941 7864487986 78e - 882, - 1 , 72989021 00 4 3226e - 882, 1 , 775 

-6. 859448788578791e-B82, 8, 97831 ' 3256 4 4 826 6e- 882, 3. 8145 

I - 6. 859448788578791e - 882, - 6, 788468549771 7 3 Be - 882, 2, 115 

--882, 2,734 ?5232e-883,l,2 7 875 



f*_ lftL _ "I 


Load P rogram [ ? J [ X J 




Look in: | fcJUJ Debug 


-1 o e &■ muh 




1^1^ 




File name: |aic23.out 


Open 




Files of type: | K .out 


■w | Cancel 




Help 
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10. Go to Debug and click on run to run the program. 

TTTTW.+JilCFirTU 

] File Edit View Project 



GEL Option Profile lools PBC DSP/BIOS Wi 



i\ & 



new linear, pjt 
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?} 



H 



Q1B3I 



§ Files 
S"Q GEL files 
B-Q Projects 
B-jE^ new lin 

Q Depend 

P~l Docurne 

SP/BIC 

General: 

Include 

f^ Libraries 

■i] rts6 

ll Source 




Breakpoints... 
Probe Points... 



Step Into 
Step Over 
Step Out 



Halt 

Animate 

Run Free 

Run to Cursor 

Set PC to Cursor 

Mu[tiple Operation 



AssemblyjSource Stepping 



Reset CPU 
Restart 
Go Main 



Reset Emulator 
Disconnect 
Restore Debug State 



ja 



Fll 
F10 
Shift+Fll 




5hift+F5 

Alt+F5 

Ctrl+F5 

Ctrl+FlO 

Ctrl+5hift+F10 



Ctrl+R 

Ctrl+5hift+F5 

Ctrl+M 



Ctrl+Shift+R 
Alt+C 



Enable Thread Level Debugging. 

Real-time Mode 

Enable Rude Real-time Mode 



Flush Pipeline on Halt 

■ UUUUMUH UUUUUb^ 
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1 1 . To see the Graph go to View and select time/frequency in the Graph and give the 
correct Start address provided in the program, Display data can be taken as per user. 
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1 1 .a) To plot the input signal: 





S3 Graph Property Dialog 


*l 


g] File Edit 


^B Project Debug t.LL Option Prohle lools VtiL D; 


^ Single Time 
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• Standard Toolbar 
GEL Toolbar 

i/ Project Toolbar 

• Edit Toolbar 
Advanced Edit Toolbar 

• Layout Toolbar 
</ Status Bar 

Debug Toolbars ► 
Plug-in Toolbars ► 




Graph Title 


Graphical Display 






new linear p|i 
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InSigBuffer 
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Disassembly 
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Plot Data From 
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Left-shifted Data Display 


Yes 




Registers ► 
Peripherals 


Autoscale 
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DCValue 







rt 


Graph ► 
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DO 

41 

JO 


Axes Display 


On 




£ 


Threads ► 


Constellation,,, 
Eye Diagram,,. 


Time Display Unit 


s 




* 


Watch Window 
Quick Watch 
Call Stack 
Expression List 
i/ Output Window 
Symbol Browser 


Status Bar Display 
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JLIlILT/l!44 UUUUU 


Magnitude Display Scale 


Linear 
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1 1 .b) To plot the noise. 
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ll.c) 



To plot the noise-corru pted signal. 
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Experiment 7: Impulse Response 
Procedure to create new Project: 

1. To create project, Go to Project and Select New. 

B^SEEBSBIBBEB] 



fflli 



File Edit View 



1 Q£ 



r 






-0 

m 
m 

IS] 



^r 



£ Files 
+ CD GE 



Debug GEL Option Profile Tools PBC DSP/BI05 V 




Use External Makefile. 
Export to Makefile,,, 
Add Files to Project,., 
Save 



Source Control 



Compile File 
Build 

Rebuild All 
Stop Build 
Build Clean 



Ctrl+F7 
F7 



Configurations,.. 
Build Options.., 
File Specific Options.,. 
Function Level Options. . , 



Project Dependencies. , . 
Show Project Dependencies 
Show File Dependencies 

Recent Project Files 



2. Give project name and click on finish. 





Project Name 
Location: 
Piujeut Type. 
Target 








r 








| C:\CCStudio_v3.1 \MyPrcjects\ 


J 


| Executable (.out) 


J 






|TMS320CG7>O< 


J 
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( Note: Location mus t be c:\CCStudio_v3.3\MyProjects ). 

3. Click on File ' ■> New ' ■> Source File, To write the Source Code. 
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AIM : To find Impulse response of a first order and second order system. 

A discrete time system performs an operation on an input signal based on predefined 
criteria to produce a modified output signal. The input signal x(n) is the system 
excitation, and y(n) is the system response. The transform operation is shown as, 



x(n) y(n) = T[(x(n)] 

► 




The convolution sum can be represented by, y(n) = x(n) * h(n) 

For Example let's find out an impulse response of a difference equation. 

The general form of difference equation is, 

M M 

y(n) = ^a k y(n-k) + ^b k x(n-k) 

k=l k=0 

Find out the impulse response of second order difference 

equation. 

Program: 

#include<stdio.h> 

#defme Order 2 

#defme Len 5 

float h[Len] = {0.0,0.0,0.0,0.0,0.0},sum; 

void main() 

{ 

int j, k; 

float a[Order+l] = {0.13 11, 0.2622, 0.13 11}; 

float b[Order+l] = {1, -0.7478, 0.2722}; 

forG=0;j<Len;j++) 

{ 

sum = 0.0; 

for(k=l; k<=Order; k++) 

{ 

if(G-k)>=0) 

sum = sum+(b[k]*h[j-k]); 

} 

if (j<= Order) 
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hU] = a[j]-sum; 
else 

h[j] = -sum; 
printf("%f ",h[j]); 

} 

} 

Output: 

0.131100 0.360237 0.364799 0.174741 0.031373 

For first order difference equation. 
Program: 

#include<stdio.h> 

#defme Order 1 

#defme Len 5 

float h[Len] = {0.0,0.0,0.0,0.0,0.0} 9 sum; 

void main() 

{ 

int j, k; 

float a[Order+l] = {0.1311, 0.2622}; 

float b[Order+l] = {1, -0.7478}; 

for(j=0;j<Len ; j++) 

{ 

sum = 0.0; 

for(k=l; k<=Order; k++) 

{ 

if((j-k)>=0) 

sum = sum+(b[k]*h[j-k]); 

} 

if(j<=Order) 

h[j] = a[j]-sum; 
else 

h[j] = -sum; 
printf("%f ",h[j]); 

} 
} 

Output; 

0.131100 0.360237 0.269385 0.201446 0.150641 
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4. Enter the source code and save the file with ".C" extension. 
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} 
if (j<=0rder) 

{ 

h[j] 
} 



a[j ]-sum; 



' 



5. Right click on source, Select add files to project. And Choose " .C " file Saved before. 
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6. Right Click on libraries and select add files to Project., and choose 
C:\CCStudio_v3.3\C6000\cgtools\lib\rts6700.1ib and click open. 
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7. a) Go to Project to Compile. 

b) Go to Project to Build. 

c) Go to Project to Rebuild All. 
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8. Go to file and load program and load ".out" file into the board.. 
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9. Go to Debug and click on run to run the program. 
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10. Observe the output in output window. 
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1 1 . To see the Graph go to View and select time/frequency in the Graph, and give the 
correct Start address or the output variable name provided in the program, Display data 
can be taken as per user. Select DSP data type as 32-bit floating point. 
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12. Green line is to choose the point, Value at the point can be seen (Highlighted by 
circle at the left corner). 
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